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ABSTRACT 


An  Integral  equation  for  the  electromagnetic  field  within  a  die¬ 
lectric  body  is  given.  The  equation  Is  set  up  for  numerical  solution 
for  the  case  of  thin-wall  cylindrical  dielectric  shells  having  finite 
length. 

The  solution  of  the  integral  equation  utilizes  a  truncated  double 

t 

Fourier  expansion  of  the  field  in  the  shell.  The  integral  equation  is 
then  enforced  at  enough  points  within  the  shell  wall  to  obtain  a  suf¬ 
ficient  system  of  linear  equations  in  the  unknown  expansion  coefficients 
of  the  field.  Numerical  integration  over  the  shell  volume  is  used  to 
obtain  the  coefficients  in  the  system  of  linear  equations.  The  sys¬ 
tem  of  equations  is  solved  numerically  for  the  expansion  coefficients 
of  the  field  in  the  shell.  Calculation  of  the  backscattered  fields  and 
the  backscattering  cross  section  are  then  performed. 

A  comparison  of  the  calculated  and  measured  backscattering  cross 
section  is  made  for  rings  with  arbitrary  plane  wave  incidence  and  for 
tubes  with  axial  plane  wave  incidence.  The  agreement  is  excellent  in 
all  cases  considered. 

The  numerical  methods,  experimental  arrangement,  computer  programs 
and  suitable  extension  of  this  work  are  discussed. 

/■- 


ill 


CONTENTS 


ACKNOWLEDGMENTS  . 

VITA  . 

LIST  OF  TABLES  . 

LIST  OF  FIGURES  . 

Chapter 

i 

I  INTRODUCTION  . 

II  FORMULATION  OF  THE  INTEGRAL  EQUATION  ... 

A.  The  Equivalent  Currents 

B.  The  Integral  Equation 

III  SOLUTION  OF  THE  INTEGRAL  EQUATION  . 

A.  The  Field  Expansion  and  the  System 

of  Linear  Equations 

B.  Incident  Field  Expansion 

IV  INTEGRATION  PROCEDURES  . 

A.  General 

B.  The  Field  at  the  Center  of  the 

Singular  Cell 

C.  The  Far  Scattered  Fields 

V  COMPARISON  OF  CALCULATED  AND  IOSURED 

RESULTS  . 

A.  Dielectric  Ring  Backscattering 

B.  Dielectric  Tube  Backscattering 

VI  CONCLUSIONS  . 

Appendix 

A  EXPERIMENTAL  METHODS  ANO  EQUIPMENT  . 


Page 

11 

iii 

v  11 

vi  11 


1 

5 

5 

7 

9 

9 

19 

23 

23 

25 

31 

36 

36 

50 

54 


57 


1v 


CHAPTER  I 


INTRODUCTION 

It  Is  a  relatively  difficult  task  to  accurately  calculate  the  scat¬ 
tering  characteristics  of  dielectric  bodies  when  they  are  placed  in  an 
arbitrary  electromagnetic  field.  Knowledge  of  dielectric  scattering  is 
important  when  a  protective  dielectric  shell  is  placed  over  an  antenna 
and  when  radar  mapping  is  used  for  geophysical  or  military  purposes. 

The  research  discussed  in  this  paper  was  motivated  primarily  by  interest 
in  the  scattering  from  protective  dielectric  shells  (radomes)  when  placed 
over  an  antenna. 

By  comparison  with  the  effort  expended  on  electromagnetic  scat¬ 
tering  from  perfectly  conducting  bodies,  the  amount  of  research  on 
scattering  from  dielectric  bodies  is  small.  Some  of  the  earliest  work 
on  dielectric  scattering  was  done  by  Lord  Rayleigh.  Stratton  [13 
points  out  that  Rayleigh  applied  the  electromagnetic  theory  of  light 
to  scattering  by  dielectric  bodies  which  are  small  in  comparison  with 
a  wavelength.  Lord  Rayleigh  used  his  results  to  describe  scattering 
by  colloidal  particles  and  to  explain  the  blue  cast  of  the  sky. 

Rigorous  solutions  for  plane  wave  scattering  from  dielectric  bodies 
have  been  found  for  the  sphere ,£2]  the  infinite  circular  cylinder, £3, 43 
the  infinite  elliptic  cylinder ,£5]  the  Infinite  parabolic  cylinder£63 
and  the  Infinite  plane  dielectric  slab .£73  Tice  and  Adney[83  for¬ 
mulated  a  rigorous  solution  for  a  dipole  within  a  spherical  shell  and 
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AndreasenC9]  developed  an  asymptotic  solution  for  a  dipole  within  a 
thin  spherical  shell. 

The  homogeneous  vector  Helmholtz  equation  must  be  satisfied  in  a 
source- free  homogeneous  medium  by  the  vector  potential  functions  for 
the  electromagnetic  field.  Morse  and  FeshbachClO!]  have  observed  that 
this  equation  is  separable  in  only  six  coordinate  systems  (rectangular, 
circular  cylindrical,  elliptic  cylindrical,  parabolic  cylindrical, 
spherical  and  conical).  They  also  point  out  that,  although  rigorous 
series  solution  can  be  obtained  for  the  remaining  five  coordinate  systems 
In  which  the  scalar  Helmholtz  equation  is  separable, Cl lH  "the  fitting  of 
boundary  conditions  Is  well-nigh  impossible  of  attainment."  The  number 
of  problems  for  which  a  rigorous  series  solution  is  practical  Is  there¬ 
fore  severely  restricted. 

Various  approximate  methods  have  been  used  to  study  the  scattering 
from  finite  dielectric  bodies.  CohenCl2D  developed  approximate  scat¬ 
tering  formulas  based  on  the  reaction  concept  and  applied  these  for¬ 
mulas  to  the  infinite  dielectric  cylinder.  Montroll  and  Hart[133 
derived  formulas  for  the  scattering  of  finite  cylinders  by  approximating 
the  fields  in  the  finite  cylinder  by  those  in  the  infinite  cylinder. 
Llnd[14]  also  used  the  infinite  cylinder  field  to  approximate  that  in 
a  finite  cylinder  but  added  a  "normal  mode"  correction  to  this  solution. 
OguchIDS]  derived  formulas  for  the  scattering  by  dielectric  spheroids 
of  small  eccentricity  by  considering  the  spheroids  to  be  perturbations 
of  a  sphere.  Ph1111pson[16]  calculated  the  scattering  from  dielectric 
rings  using  an  Iteration  method.  Geometrical  optics  was  applied  by 
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Peters  and  Thomas[17]  to  scattering  from  thin-wall  spherical  shells  and 
by  Kouyoumjian,  et.al.[18]  and  Peters,  et.al.[19]  to  more  general  die¬ 
lectric  bodies. 

For  the  Scattering  from  an  arbitrary  dielectric  body,  rigorous 
series  methods  are  impractical  and  the  accuracy  of  the  approximate 
techniques  is  unknown.  For  these  reasons,  another  approach,  the  in¬ 
tegral  equation  method,  has  been  applied  in  recent  years  to  some 
dielectric  scattering,  problems .  The  integral  equation  method  is  fun¬ 
damentally  a  rigorous  numerical  technique  for  finding  the  fields  within 
a  dielectric  volume. 

The  power  of  the  integral  equation  method  was  recognized  in  an 
early  report  by  Rhodes. [20]  A  recent  book  by  Harrington[21]  devotes 
a  considerable  amount  of  discussion  to  the  solution  of  integral  equa¬ 
tions  for  field  problems  by  means  of  the  method  of  moments.  Richmond 
has  used  the  integral  equation  method  successfully  to  solve  a  varied 
group  of  dielectric  scattering  problems:  calculation  of  radome  dif¬ 
fraction  patterns ,[22]  diffraction  by  metallic  or  dielectric  toroids 
with  a  coaxial  magnetic  line  source,[23]  scattering  from  finite  die¬ 
lectric  cyl 1nders[24]  and  scattering  from  infinite  dielectric  cylinders 
of  arbitrary  cross  section. [25, 26]  Waterman[27]  has  recently  applied  an 
"extended"  Integral  equation  method  to  scattering  by  dielectric  bodies; 
he  utilizes  spherical  mode  expansions  and  satisfies  the  integral  equa¬ 
tion  over  the  body's  interior  region. 

The  solution  of  an  integral  equation  for  the  field  within  a  body 
Is  usually  accomplished  by  approximating  the  equation  by  a  system  of 
linear  equations.  The  unknowns  in  the  system  are  either  the  field  it- 
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self  at  a  number  of  points  within  the  body  or  a  set  of  expansion  co¬ 
efficients  for  the  field  in  the  body.  Solution  of  large  systems  of 
equations  and  repeated  computation  of  complicated  volume  integrals 
require  large,  high  speed  computers.  It  is  for  this  reason  that  the 
integral  equation  method  was  not  applied  extensively  prior  to  the  ad¬ 
vent  of  the  modern  generations  of  digital  computers. 

In  general,  a  radome  can  be  an  arbitrary  three-dimensional  shape/ 
may  be  inhomogenepus  and  Is  placed  in  the  relatively  complex  near  field 
of  a  radar  antenna.  Solution  of  such  a  difficult  scattering  problem 
must  proceed  in  a  sequence  of  smaller,  yet  significant,  steps.  Thus,  as 
the  first  step  in  this  paper,  the  integral  equation  method  is  applied  to 
the  problem  of  calculating  the  scattering  by  thin-wall  dielectric  cir¬ 
cular  cylindrical  tubes  of  finite  length.  In  Chapter  II  the  equivalent 
source  currents  for  the  scattered  field  and  the  basic  integral  equation 
are  discussed.  The  expansion  of  the  field  in  the  shell,  the  derivation 
of  the  system  of  linear  equations  and  the  Fourier  expansion  of  an  in¬ 
cident  plane  wave  are  described  in  Chapter  III.  The  techniques  of  nu¬ 
merical  integration,  the  special  manner  of  integrating  through  the  singu¬ 
larity  and  the  far  field  calculation  are  discussed  in  Chapter  IV.  In 
Chapter  V  a  comparison  is  made  of  the  calculated  and  measured  plane-wave 
backs cattering  from  dielectric  rings  and  cylindrical  shells.  Conclusions 
are  presented  In  Chapter  VI.  A  discussion  of  the  experimental  method 
and  a  description  of  the  computer  programs  used  are  presented  in  the 
Appendices. 


CHAPTER  II 


FORMULATION  OF  THE  INTEGRAL  EQUATION 
A.  The  Equivalent  Currents 

The  following  equivalence  principle  and  its  derivation  were  de¬ 
veloped  by  Rhodes. [203  The  derivation  below  is  the  same  as  that 
given  by  Richmond. [283 

The  time  dependents  e+^wt  is  understood  and  linear,  non-magnetic 
(y*u0)»  Isotropic  dielectric  media  are  assumed.  The  medium  may  be 
inhomogeneous  and  lossy,  i.e.,  it  may  have  a  complex  permittivity. 

Let  a  current  source  J(x,y,z)  in  a  medium  (y^x.y.z),  e^x.y.z) 
generate  a  field  (E^ (x»y ,z) ,  H^x.y.z).  This  field  will  satisfy 
Maxwell's  curl  equations: 

(1)  x  Hj  *  J  +  j  we,j  Ej 

)  v.  x  E.  =  - j  u>.  ui  Hj 

let  the  same  current  source  J (x,y,z)  in  a  new  medium  (u  (x,y,z),  e(x,y,z)) 
generate  a  different  field  (E.(x,y,z),  H(x,y,z)).  This  new  field  must 
also  satisfy  Maxwell's  curl  equation: 


(2) 


7  x  H  8  J  +  j  id  e  E 
7  x  E  =  -j  u  g  H 
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The  difference  between  these  fields  is  called  the  scattered  field: 


If  Maxwell's  equations  for  the  two  sets  of  fields  are  subtracted,  we 

obtain  the  curl  equations  for  the  scattered  field: 

! 

(4)  fv  x  =  j  w(e  -  e.)  E  +  j  to  ei 

X  Eg  =  -j  u>(u  -  ^0  H  ‘  3  Hj 

The  term  j  w(e  -e.)  E.  (which  is  nonzero  only  where  e  f  e^)  may  be 
Interpreted  as  an  equivalent  electric  current  density  radiating  in  the 
original  medium,  (y^  e^).  Similarly,  the  term  -j  cu(y  -  y.)FL,  can  be 
Interpreted  as  an  equivalent  magnetic  current  density  radiating  in  the 
original  medium  (yj.e,.).  These  current  densities  are  the  sole  source: 
of  the  scattered  field. 

In  this  investigation,  the  original  medium  is  chosen  to  be  free 
space  so  *  eQ.  Furthermore,  all  media  are  assisned  nonmagnetic  so 
that  Uj  *  u =  uQ  and  the  equivalent  magnetic  current  density  vanishes. 
The  electric  current  density,  j  w(e  -  eQ)  E^,  then  becomes  the  sole 
source  of  the  scattered  field  and  radiates  in  unbounded  free  space. 

As  noted  earlier,  the  equivalent  current  concept  can  be  applied 
to  inhomogeneous,  lossy  dielectric  media. 
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B.  The  Integral  Equation 

The  fields  radiated  by  a  source  current  density  0  in  free  space 
can  be  calculated  using  an  integral  form  such  as 

(5)  Ep  =  ///  J  •  Gp  (r,r')  dv’  , 

where  Ep  Is  the  p-th  component  of  the  radiated  field,  is  the  source 
current  density,  V'  is  the  source  volume  and  represents  the  free 
space  vector  Green's  function  for  the  p-th  component.  The  equations 
given  by  Richmond^ 8D  and  given  in  Eq.  (9)  can  readily  be  put  into  the 
above  form.  Substituting  the  equivalent  current  density  of  Eq.  (4) 
Into  Eq.  (5)  and  using  this  result  in  Eq.  (3),  we  obtain  our  integral 
equation  for  the  field  in  the  dielectric  body  and  the  exterior  free- 
space  region: 

(6)  EP  =  Ejp  +  JJj  j  «(e  -  e0)  E  •  ^  (r,r')  dv' 

To  determine  the  field,  it  is  sufficient  to  enforce  Eq.  (6)  in  the 
interior  region. 

Equation  (6)  must  be  satisfied  for  each  of  the  three  components 
of  the  field  such  as  Ex,  Ey,  and  Ez  in  rectangular  coordinates 
or  Ep,  E*  and  Ez  in  the  cylindrical  system. 

Solution  of  Eq.  (6)  yields  a  solution  for  the  total  field  within 
the  body.  Then  the  field  at  any  exterior  point  can  be  calculated  from 
Eq.  (6).  For  the  scattering  by  cylindrical  dielectric  shells,  a  tech- 
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,  •  in«v  and  lossless,  homogeneous  an 

geneous  ct  our  attention  to  a  specific 

trate  the  techniques  we  ere  ,ar  cylindrical  shells. 

case:  thin-wall .  homogeneous ,  dele  ^  any 

The  shell  is  assumed  homogeneous  and  .uffic 

.  „f  the  fields  is  negligible  within  the  shell, 
radial  variation  o  c.  ,  At  any  point, 

Kell  and  its  coordinate  system  are  illustrated  ,n  Fig.  1-  » 

I  „  the  field  components  can  be  expanded  in  a  Fourier 
*.  *’""9  the  shell,  the  field  c  P  are  functl0„5 

in  a  The  expansion  coefficients  i 

series  in  <t>.  Tne  e  v  \  Both  series 

of  t  and  are  likewise  expanded  in  a  Fourier  senes  0 

are  truncated  to  finite  su.  as  a  practical  matter. 

truncation  depends  on  the  particular  shell  and  the  me  ■ 

,  ts  the  shell  to  obtain  a  set  of  simultaneous  linear 
at  many  points  in  ^  ^  so)>ed  for  the 

equations  in  the  expansion  «■»"«’"  •  ^  ^  wUM„  the  shell. 

expansion  coefficients,  or  e,u„a  ,  ^  awtotlc 

a.  ah.  iTitai  field  vnthm  tne  snei  > 

Finally,  when  the  total  e 

form  of  Eg.  (6)  is  employed  to  calculate  the  dis 


CHAPTER  III 

SOLUTION  OF  THE  INTEGRAL  EQUATION 
A.  Field  Expansion  and  the  System  of  Linear  Equations 

The  coordinate  system  and  a  cylindrical  shell  are  shown  in  Fig.  1. 
The  circumflex  Indicates  a  unit  vector. 


Fig.  1—  Cylindrical  coordinate  system 
and  cylindrical  shell. 

The  shell  Is  assumed  to  be  sufficiently  thin  that  field  variations 
with  p  can  be  neglected  within  the  shell.  In  the  dielectric  region, 
the  fields  can  be  expanded  in  a  Fourier  Series  as  follows: 
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(7) 


f  oo 

Ep(z)  *  l  Ift  (z)  cos(n<j>)  +  (z)  sin  n$] 

n=0  n  n 

<  E*(z)  *  l  CAj  (z)  cos(n*)  +  B*  (z)  sin  n*3 
n=0  M 

EZ(z)  s  I  CA  Z(z)  cos  n4>  +  B*  (z)  sin  n<j>3 
n=0  " 


In  Eq.  (7),  as  noted,  the  An's  and  Bn's  are  functions  of  position 

along  the  shell  and  can  be  represented  by  a  Fourier  series  in  z.  The 

period  of  the  expansion  is  p  where  p  is  greater  than  or  equal  to  the 

shell  length.  The  expansions  of  the  A  's  and  B  's  follows: 

n  n 


(8) 


■S 


(z)  *  Jo  Cb“”"  “!  ^ z)  +  bs”" 5i"  (2r 2)3 
<2>- Jo^  »  (f-zl  +  asl. stn  (?r 2)3 

(z>  =  Jo cb™ cos  (2r z)  +  bs”" sin  <2r 2,3 


A"Z(Z)  =  JoCa“"  “s  <2fz>t<s1n 

B"  (Z)  '  ml  Cb*"  C°S  Z>  +  bS""  Sl"  ‘T"  2)3 

km 

In  the  circular  cylindrical  system,  there  will  be  no  coupling  be¬ 
tween  modes  of  different  index  n  in  the  ^-expansion.  For  this  reason 
It  is  possible  to  solve  the  Integral  equation  for  one  value  of  n  at  a 
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time,  to  obtain  the  solution  for  all  necessary  n  and  then  to  superimpose 
the  results  to  determine  the  total  field  within  the  shell.  This  allows 
a  considerable  saving  In  computer  storage  compared  with  solving  the 
entire  system  in  one  fell  swoop.  The  z-expansion  must  be  considered 
In  Its  (truncated)  entirety  for  each  value  of  n,  however. 

Eq.  (6)  is  to  be  enforced  at  a  number  of  points  in  the  shell.  For 
numerical  purposes,  the  shell  is  divided  into  equal-length  elemental 
rings,  each  ring  being  short  compared  with  the  wavelength.  The  numerical 
Integrations  over  the  tube  (shell)  are  performed  on  a  ring-by-ring  basis 
and  the  total  Integral  is  found  by  summing  the  contributions  of  all  of 
the  rings.  Fig.  2  illustrates  the  subdivision  of  the  cylindrical  shell. 


z 


Fig.  2—  Division  of  shell  into  elemental  rings. 
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The  match  points  (points  of  enforcement  of  the  integral  equation) 

are  at  the  center  of  the  wall  cross  section  of  the  matching  rings  (two 

per  ring)  and  are  located  at  <f>= 0  and  ♦*!  -.  The  matching  rings  are 

n  2 

evenly  spaced  along  the  tube.  If  M  is  the  maximum  index  of  the  z-expan- 
sion,  (2M  +  1)  matching  rings  are  required  to  obtain  a  complete  system 
of  equations  in  the  field  expansion  coefficients. 

The  next  step  is  to  evaluate  the  integrals  for  the  scattered 
field  at  points  in  the  shell.  In  the  following  work,  attention  is 
confined  to  the^n-th  mode  in  $*,  it  is  understood  that  in  general  many 
such  modes  may  be  needed  to  obtain  an  accurate  solution  for  the  total 
field  in  the  shell. 


The  following  equations  are  given  by  Ri chmondE 2PQ  for  the  cartesian 
components  of  the  field  generated  in  free  space  by  a  source  current 
density  with  components  Jx,  Jy,  and  Jz: 


(9) 


JJJ  e"J’kr  (P(r)  Jx  +  (x-x')Q(r)C(x-x’)Jx 
v1 

+  (y-y')Jy  +  (z-z')  JZ]}  dv' 

fff  e'jkr  {P(r)  Jy  +  (y-y’)Q(r)[;(x-x,)Jx 
v' 

+  (y-y')Jy  +  (z-z’)Jz3}  dv* 

Iff  e’jkr  {P(r)  0Z  +  (z-z')Q(r)[(x-x')Jx 
v1 

+  (y-y')Oy  +  (z-z')J2]}  dv' 
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In  Eq.  (9), 


In  Eqs.  (9)  and  (10),  k  =  o>  ^oeo  -  2tt/a  and 

r  =  p2  +  p'2  -  2pp*  cos(  <M‘)  +  (z-z1)2 
is  the  distance  from  the  source  to  the  observation  point.  The  equiv¬ 
alent  current  density,  joj(e-e0)£,  is  substituted  into  Eq.  (9)  to  ob¬ 
tain  an  expression  for  the  scattered  field  in  terms  of  the  total  field. 
The  field  components  are  then  converted  to  circular  cylindrical  com¬ 
ponents  and  the  field  expansion  of  Eq.  (7)  and  (8)  used  to  get  the 
following  result  for  the  n-th  mode  scattered  field  where  cm  *  cos(-?!5l  2), 

III  P 

sm  *  sin(~r*)  and  er  -  e/cQ. 


I  ["{(a’  +  as|L  sjcos  n$’  +  (b’_  c+bs’sm)sin  n<j> ’  }cos<j>' 

<P(r)  mn  m'  mn  m  mn  m 

}  -{(a*  c+as*  sj  cos  n$'  +  (b^c+bs*  sjsin  n4>’  >sin<j>  ' 

I  v  mn  m  mn  m  T  mn  m  mn  m 


+Q(r)(x-x')2 


{(aJLc L+asJLs)  cos  n$'  +  (b*Lcm  +  bs^  sjsin  n^' >  cos*1' 
nmn  m  mn  m  mn  m  mn  m 


{(a*  cm  +  as’  s_)  cos  n$'  +  (b’  cm  +  bs’  s  )sin  n$' }sin$‘ 
. 'Tnn  m  mn  m  r  mn  in  mn  nr 
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+  Q( r) (x-x ' ) (y-y  ‘ )  f  {(a®ncm  +  as®  s_)  cos  n*'  +  (b®  c  +  bs®  s_) 
^  ^  mn  m  mn  m  mn  mn  mn  m 

sin  n$' }  sin$‘ 

(11)  +f(a*  (L  +  as*  s  )  cos  n<*> •  +  (b*  c  +  bs*  s_) 

(cont  )  a-  ran  in  mn  m  ran  ra  mn  m 


sin  n$* }  cos  #' 


+  OfrMx-x'itJ-i1)  +  as^nsm)  cos  m' 

*  (bmoc«,  *•»»*•>  s1"  "♦] 


c  y  _  cr  ~1 

sn  “  4n 


L  n!  e 

m=0  v 


{ ( a^n c  +  as®  sm)  cos  n*'  +  (b®  c_  +  bs®  sjsin  n«i>'}sin  4> * 
mn  ra  ran  ra  ran  ra  ran  m  r 


sin  n<(> ' }  cos*' 

L^amnSjl  +  “mnsm>  cos  *  ^nS.  +  bsmnsm> 

sin  n<j.' }  sin  <*»' 

9  T  ( ( a*  c  +  as®  sm)  cos  n<t>  ’  +  (b®  c_  +  bs®  sj 
Q(r)  (y-y1)2  mn  ra  ran  ra  ran  ra  ran  ra 

sin  n<j>'}  sin  4> ' 

+{(a*  cm  +  as*  s_)  cos  n$'  +  (b*  c  +  bs*  s  ) 

__  mn  ra  ran  ra  mn  m  mn  m  j 

sin  n<fr'}  cos<j>' 


+t(aL«cm  +  asm«sm)  cos  n^‘  +  ( b *  c  +  bs*  sj  sin  n^*  >cos<j> ' 
mn  m  ran  ra  mn  m  mn  ra 


{(a®c  +  as®  s_)  cos  n<j>'  +  (b®  cm  +  bs®  s  ) 
+  Q(r)(y-y')(x-x')  mn  m  mn  m  mn  m  mn  m 


*  Q  (r)(y-y')(z-z')  <V*  +  aW  cos  »♦'  +  (b,^  *  bs^sj  I  dv' 

L  sin  n*’  J  J 
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e-  -1 


(13)  ESJ  •  -fc-  I  ///  e*ikr. 

sn  ’  m>0  v' 


P(r)  C<4cm  +  asmnsm>  cos  +  <bmnS*  +  bsmnsm>  si"  "*'3  . 


+  Q(r)(z-z')(x-x') 


Ha  *  <a)u* *Ka  *  b%L*J 

sin  n<t»' }  cos  $’ 

^{(a^nC,-  +  as*  s  )  cos  n$'  +  (b*  cm  +  bs*  s_)  J 
u  mn  m  mn  m  mn  m  mn  in' J 

sf  n  n  }  sin 


+  Q(r)(z-z')(y-y') 


r{(a£n(L  +  asL,sm)  cos  n<j>'  +  (b£  c  +  bsfLsJn 
mn  m  mn  m  mn  m  mn  m  I 

sin  n$' }  sin  <t>' 

±^alncm  +  aslsJ  cos  n$'  +  (b*  c  +  bsJLs  )  . 
“  mn  m  mn  m  mn  m  mn  m' J 

sin  n  $' }  cos 

.Z  „  A  u_z 


*  Q(W(z-2’)2  [(«i,Si  +  cos  *  (ba„cm  +  bsansro)  sin  n*'|  }dv' 


■]} 


The  symmetries  of  the  integrands  for  the  scattered  fields  about 

*  3  0  and  #  *  jj- 1-  permit  a  convenient  correspondence  between  the  inte¬ 
grations  for  a  match  point  at  4>  *  0  and  one  at  <t>  3  Only  the  inte¬ 
gration  for  <p  *  0  need  be  calculated  and  the  ones  for  <j>  =  ^  \  can  be 
obtained  from  them.  Also,  the  odd  symmetry  of  some  of  the  integrand 
terms  results  in  those  terms  integrating  to  zero  for  a  match  point  at 

♦  *  0  or  ♦  s  ^  |  •  After  making  use  of  the  odd  symmetry  of  certain 
terms  and  of  the  relationship  between  the  $  *  0  and  <j>  ■  jj-  j  calculations, 
the  cylindrical  components  of  the  scattered  field  for  a  given  matching 
ring  can  be  written: 
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Es>-*>  ■  <v»  i. 

m=0 


f  3L?n  Al  (z)  , 
mn  mn  + 


af  A2  (z)  +  b  A3  (2) 
mn  mn'  7  mn  mn 


asL  B1 _ ^  +  asL  +  bs*  B3  (z)| 


mn  mn 


mn  mn 


mn  mn 


(l^E^o.z)  =(Ey.-l)  l 


s  J 

b°  A4  +  bz  A5  +  a1^  Afi  ( 7 1  1 

mn  mn  Dmn  ASmn  A6imrz' 


m=0 1 


+  bsmn  B4  ^  +  bs2  B5  ^  +  as^  B6  (z) 
mn  mn  mn  '■'mn  u°m.>  '•  > 


mn  mn 


mn  mn 


“mn  AV<Z)  *  A8-(Z)  +  bl  A9„(z) 


mn  mn 


E„(°.i)  •  (yi)  I  J  mn  mn 

m r  asmn  +  as2  SB  ^  +  bs ^  B9  (z 

1  mn  mn  mn  mn  mn  mn 


mn  mn 

>2  B8 
mn  mn 


r 


(15X 


n  1  ir  00  f  bp  Al  +  b4  A2  ' z ^  -  a'*1  A"?  fz) 

Esn  (nl’z)=  <Vl)  I  I  ""  mn  m"  ™  J 

"■“O’-  -  ,z)  '  '  2  <z>  -  asi  B3^z| 


bs  B2 
mn  mn 


mn  mr 


Esn  rz>"  <V‘>  I  1 

m=0 


(?)  +  (jz  A9  (Z) 

\  Kn  «„(2)  * 

f-»l  «J2)  -  4n  A5»n<2)  +  «JZ> 

(2>  '  «*  B5J2)  *  bsm„  86miz: 


Esa  <H’Z>*  «V«  |  J  "'mn 
mV 


-as  B4 
mn  mn 


b!L  A7 _ ^  +  b2  aa  (z) 


'mn  A8mn 


+  bs£  B7mn(z)  +  bs2  __ 
mn  mn  mn  mn 


-  a/'  A9  iz) 
mn  mrt 

bs‘  B8  (z) 


as?  B9  (z 
mn  mrt 


Almn(z)  through  A9mn(z)  and  B1fnn(z)  through  B9^n ( z )  are  defined 
below.  The  Integrals  are  to  be  calculated  for  an  observation  point 
at  the  center  of  the  shell  wall  cross  section  for  $=  0  and  z  equal  to 
the  z-coordinate  of  the  match  point  in  question. 
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A1mnZ^  "  4ff  //f  e”^kr  <p(r)  cos  ♦'  +  Q(r)(x-x')2  cos  *' 

+  Q  (r)(x-x')(y-y') }  cos  z')  cos  n$‘  dv' 

A2mnZ)  =  ~  ///  e”jkr  {Q(r}(x-x‘)(z-2')  }  cos  (^V)  cos  n*'dv' 

4tt  v'  v 

A3mJZ>  *  b  //f  e  "jkr  <-*»*)  s1n  *'  "  Q^Kx-x')2  sin  ♦' 

+  Q(r)(x-x')(y-y ')  cos  <j>'}  cos  (y  z')  sin  n$'  dv' 

A4mnZ^  a  37  J/f  e"jkr  {p(r)  s1n  ♦  '  +  Q( r) C x- x ' ) (y-y * )  cos 
+  Q(r)  (y-y')2  sin  <*.' }  cos  (^p-z')  sin  n<f>'  dv' 

A5mnZ^  “37  J/j  e  "Jkr  {Q(r)(y-y')(2”Z’)>  cos  (^p-z')  sin  n-t'dv' 

< 

) 

^mr\^  ~  37  e  ~jkr{p(r)  cos  $'  -  Q(r)(x-x')(y-y ')  sin  <f>' 

+  Q  (r)(y-y')2  cos  4'}  cos  (^p-z')  cos  n$$>  dv' 

A7m^  ■  ^  ///  e“jkr  (Q(r)  (z-z’)(x-x’)  cos  *• 

+  Q(r)(z-z‘)(y-y')  sin  41'}  cos  (^—z * )  cos  n4>*  dv* 

A8mnZ)  a  37  UJ  e‘jkr  {P(r)  +  Q(r)(z-z')2}  cos  (^z')cos  n*'  dv' 

A9mn2^  "  b  //{  e"jkr  {-Q(r)(2"2')(x-x')  sin  4' 

+  Q(r)(z-z')(y-y')  cos  4'}  cos  (^jp-z')  sin  n<t>'  dv' 
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B1mn^  through  B9mn^  are  found  directly  by  replacing  cos  (~pz') 
by  sin  z ' )  In  the  corresponding  equations  for  Al  through  A9  ^2 \ 
The  Integrations  for  the  Bmn'2  *  are  also  evaluated  at  <t>  *  0. 

The  field  expansion  of  Eq.  (8)  and  the  results  given  In  Eqs.  (14) 
and  (15)  are  used  In  Eq.  (3)  to  obtain  the  set  of  linear  equations 
corresponding  to  one  point  along  the  shell,  l.e.,  for  one  matching  ring. 
The  summation  over  m  has  been  truncated  to  a  maximum  Index  value  of 
m  »  M.  For  this  case  a  total  of  (2M  +  1)  sets  of  the  following  equa¬ 
tions  will  be  necessary  to  have  a  sufficient  system.  In  Eq.  (17),  the 
dependence  of  Almn,  etc.  upon  z  has  not  been  shown  explicitly;  It  Is 
understood  that  these  coefficients  are  functions  of  the  z-coordlnates  of 
the  matching  rings. 


•  ap  (z)+  (e  -1)  V  amnA1mn+asmnB1mn+amnA2mn 

^  "  H+asmnB2mn+bmV3mn+bslB3mn; 


(er-l)  l  b";M™>+bsmnB4»,n+bmnA5nm 


2  M  ^„A7mn+asP,B7mB+a2  A8L 

az  (z)+  (e  -1)  y  mn  ^  mn  'nn  'nr 

1n  r  m-0  Uas2  B8m +blA9mM+bsiB9 

mn  mn  mn  mn  mn  n 


b?  (z)  +  (c  -i)J  i  b|”"*1™>+ll,™B1"»>+bnmA2mii 
r  ■"-»  1  +bsmnB2mn-^nA3mn-as^B3n, 
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r 


M 

l 

m=0 


(17)< 

cont 


M 

l 

m=0 


1 *<^2jr4 


=  bin(z)  +  (V1^ 


^nA4, 


as  B4  -a  A5 
mn  mn  mn  mn  mn 


m=0  -  as  *;  B  5  +b A6  +bsi  B6 
mn  mn  mn  mn  mn  mrt 


K„si"<F*> 


M 


2m*  J  =  bin^z^+  (£r_1^  .  1 


f  b* nA7  +bs£nB7  +b*A8  1 

.»  mn  mn  mn  mn  mn  mn  I 

.2 


mn  mn  mn  mn 


m=0  ,+bs";  B8  -a^  A9  -as^  B9_  , 
[_  mn  mn  mn  mn  mn  mnl 


The  a£n(z),  a^z),  etc.  are  the  n-th  mode  expansion  coefficients 
in  the  Fourier  expansion  (in  $)  of  the  incident  field  at  the  particular 
matching  ring  in  question. 


R  Incident  Field  Expansion 

In  general,  the  incident  field  may  be  very  complex,  e.g.,  the  near 
field  of  an  antenna.  To  illustrate  the  technique,  however,  we  shall  con¬ 
sider  the  case  for  which  the  incident  field  is  a  linearly  polarized 
plane  wave  with  arbitrary  propagation  angle  e.  For  this  case  a  formula 
for  the  Fourier  expansion  of  the  incident  field  can  be  found  quite 
readi ly .  . 

The  Fourier  expansion  of  the  incident  plane  wave  over  a  ring 
of  radius  a  can  be  found  from  the  general  form  of  the  plane  wave  at  any 
point  in  space.  Fig.  3  illustrates  the  geometry. 

The  phase  of  the  incident  wave  is  relative  to  the  origin  and  the 
rectangular  field  components  have  been  converted  to  circular  cylindrical 
components  in  the  equations  below. 

For  the  TE  case,  the  incident  electric  field  is  given  by: 


TE 

LjfPfi’tZ)  =  (p  cos  $  -  $  sin  <f>)  e 


+jkpsin  <j>  sine  e+jkzcose 


(18) 


tall  XI 
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X  OBSERVATION 


Fig.  3— Linearly  polarized  Incident  plane  wave. 


(19) 


The  Fourier- Bessel  expansion  is  needed: 

QD 


Jupsin*  _  ejupcosU-  j)  =  \  N  jn  J  (up)  cos  n(4>-  5-) 
'•  c  n=0 

In  Eq.  (19),  u  =  k  sine 

"n  2  in  >  1 

0  is  the  n-th  order  Bessel  function  of  the 
n 

first  kind. 


Equation  (19)  is  used  In  Eq.  (18)  and  after  simplification,  the 
TE  Incident  field  can  be  expressed  by  the  following  Fourier  expansion: 


C^.jCuf))  +  Jn+1(up)D  cos  n 


(20)  E,(p,*,z)  =  p  (  l 

[  n=l ,3*5. . . 

+  I  j  D  ,(up)  +  0  .(up)D  sin  n$  e^k2COse 

n»2,4,6...  n  1  n+1 

+  4>j-j  Jj(up) 

+  c  'Pn-l{up)  "  Vl(upn  s1n  n* 

+,  3  s  -  Jnrt(up)3  cos  o»bkzcos9 

n=2,4,6. . .  J 

For  the  TM  case,  the  Incident  electric  field  Is  given  by: 


(21)  £j  M  *  (p  cose  si n<p  +  <p  cose  cos<t>  -  z  sine)  e 


jups1n$ejkzcose 


Comparison  with  Eq.  (18)  shows  that  with  the  exception  of  the 
2-component,  the  TM  expansion  can  be  obtained  directly  by  using  the  re¬ 
sults  of  the  TE  expansion.  The  final  result  is: 

TM  f 

(22)  £..j(p ,4>,z)  =  p  cos  -j  j  Oj(up)  +  l  CJn_j(up)-0n+1(up)Us1n  n$ 

l  n-1,3,5. . . 

+  n=2^4,6...  ’J  CJn-l(up)  *  Jn+l(up^cos  "*]  ejkzcose 

+  i  cose  (  l  C0n., (up)  +  Jn+7(uP)3  cos  n$ 
[n-1,3,5...  n  L  n+1 

^  n-2^4,6  .  j  CJ"-l(Up)  +  WUp):i  sin  n*  ejkzcose 
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(22)  -  z  sine  1 0  (ua)  +  £  20  (up)  cos  n$ 

cont.  L  0  n=2,4,6...  n 


+  I  j2  J  (up)  sin  n^le^20056 
n-1,3,5...  n  J 

The  Fourier  expansions  over  a  ring  of  radius  a  can  be  found  by 
substituting  a  for  p  in  Eqs.  (20)  and  (22). 


CHAPTER  IV 


INTEGRATION  PROCEDURES 

This  Chapter  discusses  the  techniques  of  numerical  integration  used 
In  the  solution  of  the  cylindrical  shell  scattering  problem.  The  gen¬ 
eral  method  of  calculating  the  scattered  field  in  the  shell,  the  for¬ 
mula  for  handling  integration  through  the  singular  points  and  the  for¬ 
mulas  for  calculation  of  the  far  scattered  fields  are  discussed. 

A..  General 

The  coefficients  defined  in  Eq.  (16)  must  be  evaluated  at  $  =  0 
for  each  value  of  the  4>  mode  index  n  and  the  z-mode  index  m  at  each  of 
the  matching  rings.  'The  general  procedure  for  the  numerical  integra¬ 
tion  is  to  divide  the  shell  into  elemental  rings  as  shown  in  Fig.  2. 

The  integrals  are  then  calculated  over  each  ring  for  each  of  the  match 
points  at  $  =  0  and  the  total  integration  is  obtained  by  adding  the  con¬ 
tribution  of  all  the  rings. 

Integration  over  a  given  elemental  ring  is  accomplished  by  dividing 
the  ring  into  a  number  of  subcells  as  shown  in  Fig.  4.  The  maximum 
arc-length  of  any  subcell  is  much  less  than  a  free  space  wavelength  and 
the  angle  subtended  by  the  subcell  at  the  ring  center  is  small  compared 
with  the  smallest  period  needed  in  the  ^-expansion. 

For  the  case  of  a  match  point  far  from  a  given  ring,  the  trapezoid¬ 
al  rule  of  numerical  integration  is  applied  over  the  entire  ring  with 
the  center  of  each  subcell  providing  one  data  point.  For  the  case  of 
a  match  point  close  to  the  ring,  finer  subdivision  of  the  cells  closest 
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Fig.  4— Ring  subdivision  for  numerical  Integration. 

to  the  match  point  Is  required  to  obtain  accurate  results. 

The  case  In  which  the  match  point  Is  within  the  ring  being  Inte¬ 
grated  over  Is  a  rather  special  one.  Examination  of  P(r)  and  Q(r)  In 
Eq.  (10)  reveals  that  one  can  expect  severe  numerical  computation  prob¬ 
lems  as  r  approaches  zero.  In  this  case,  the  usual  numerical  Integration 
procedure  (Including  finer  subdivision  for  cells  near  the  match  point) 
is  used  to  calculate  the  contribution  of  all  the  parts  of  the  ring  ex¬ 
cept  for  that  of  the  cross-hatched  cell  In  Fig.  4;  this  cell  Is  the 
so-called  singular  cell. 

The  contribution  of  the  singular  cell  at  Its  center  Is  then  cal¬ 
culated  from  analytical  formulas.  These  formulas  are  derived  In  the 
next  section. 


B.  The  Field  at  the  Center  of  the  Singular  Cell 

Three  basic  assumptions  are  used  to  obtain  the  singular  cell  contri¬ 
bution  at  its  own  center.  First,  it  is  assumed  that  the  cell  can  be 
closely  approximated  by  a  rectangular  parallelepiped.  Second,  the  cell 
is  assumed  sufficiently  small  that  any  field  variations  within  its 
volume  can  be  neglected.  Finally,  the  linear  dimensions  of  the  singular 
cell  are  assumed  to  be  much  smaller  than  the  wavelength. 

The  rectangular  parallelepiped  is  a  suitable  approximation  to  the 
shape  of  the  singular  cell  in  an  elemental  ring  of  a  thin  cylindrical 
shel  1 . 

The  second  assumption  is  also  valid  since  the  odd  symmetry  of 
the  sin  n^  mode  about  the  center  plane  of  the  singular  cell  can  be 
shown  to  result  in  zero  field  at  the  cell  center  due  to  sin  n4>  current 
densities.  This  leaves  only  the  cos  ni>  variation  with  $  and  this  term 
has  a  zero  slope  at  the  center  plane  of  the  cell.  The  period  of  the 
z-expansion  is  assumed  large  enough  that  variations  with  z  along  the 
singular  cell  can  be  safely  neglected. 

The  singular  cell  problem  is  therefore  reduced  to  finding  the 
electric  field  at  the  center  of  a  rectangular  box  of  uniform  current 
density  flowing  in  an  arbitrary  direction.  Consider  the  box  of  cur¬ 
rent  density  to  be  aligned  with  a  local  cartesian  coordinate  system  with 
the  origin  at  the  center  of  the  box  as  shown  in  Fig.  5. 

Examination  of  the  symmetry  of  Eq.  (9)  for  this  case  shows  that 
only  Jx  will  contribute  to  Ex,  to  E^  and  Jz  to  Ez.  Thus,  the  sol¬ 
ution  needs  to  be  developed  for  only  one  component  of  current  and  a 
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Fig.  5--Rectangular  parallelepiped  of  current  density  0. 


cyclic  rotation  of  the  coordinates  will  give  the  solutions  for  the 
other  components.  The  solution  will  be  found  for  £  =  x  Jx. 

Consider  Maxwell’s  curl  equation  relating  the  electric  field  to 
the  curl  of  the  magnetic  field: 

(23)  E  •  jJ-  p  x  H  -  33 

Is  assumed  known  and  uniform  in  the  small  box.  In  order  to 
evaluate  the  field  of  this  source  at  the  center  of  the  box,  then,  v  x  H 
must  be  found  at  the  center. 

By  symmetry,  for  this  particular  problem,  there  will  be  no  magnetic 
field  at  the  center  of  the  box,  nor  wll ;  there  be  any  magnetic  field 
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whatever  at  any  point  on  the  x-axis  (for  an  x-directed  current  density). 
Thus,  for  this  case,  H  (x,o,o)  =  0  and  — ^  =  — ^  =  0.  The  curl  operation 
on  H  can  therefore  be  written  as  follows: 


(24) 


v  x 


-J 

(o,0, 


o) 


x  y  z 

0  -1  i- 

ay  az 
o  Hy  Hz 


(o.o.o) 


The  partial  derivatives  of  Hz  and  Hy  will  be  found  by  evaluating 
the  magnetic  field  at  the  center  of  the  box  (known  by  symmetry  to  be 
zero)  and  at  a  small  distances  &y  and  az  along  the  y-and  z-axes.  The 
partial  derivatives  are  then  found  from  the  definition: 


^(o.o, az)  -  Hy(o,o,o) 
AZ 


H2(o,Ay,o)  -  Hz(o,o,o) 
ay 

down  the  x-axis  at  the  box  of  current, 
can  be  seen  that  Region  I  will  not 
at  (o,Ay,o),  the  magnetic  field  at 
that  point  will  be  strictly  due  to  Region  II.  Fig.  7  shows  the  geometry 
which  is  now  considered.  For  the  special  case  when  the  distance  and  di¬ 
mensions  are  small  compared  to  the  wavelength,  the  x-component  of  the 
magnetic  field  at  a  point  a  distance  y  2  b  from  the  center  of  the  slab 
of  current  can  be  shown  to  be: 
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Fig.  6  shows  a  view  looking 
From  synmetry  considerations,  it 
contribute  to  the  magnetic  field 
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Similarly,  the  y-component  of  the  field  a  distance  2=c  along  the  z-axis 
can  be  shown  to  be: 


Reference  to  Fig.  6  reveals  that  the  incremental  fields  sought 
are  those  fields  just  calculated.  Namely, 


(28) 
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These  are  used  in  Eq.  (25)  to  find: 


Use  of  Eq.  (29)  In  Eq . (23)  gives  the  final  result  for  the  electric 
field  at  the  center  of  a  box  of  x-dlrected  current  density  Jx: 
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(30)  EX  *  A-  2J 
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Similarly,  for  Ey  and  E2, 
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If  Jx,  Jy  and  Jz  are  replaced  by  the  equivalent  currents  of  Eq. 

(4),  we  have  directly  the  scattered  fields  at  the  center  of  the  singular 
cell  contributed  by  the  current  in  that  cell: 

f  r 
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It  was  pointed  out  earlier  that  within  the  singular  cell,  only  the 
equivalent  current  component  in  the  direction  of  one  of  the  local  coor¬ 
dinates  will  contribute  to  that  component  of  the  scattered  field.  The 
above  formulas  will  therefore  be  applied  only  In  calculating  the  values 
of  coefficients  In  Eq$«(14)  and  (15)  which  represent  self-coupling  terms. 
For  example,  there  will  be  singular  cell  contributions  necessary  In  the 
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calculation  of  Almn(z)  and  81mn(z)  in  Eq.  (14)  since  these  represent 
the  "self-coupling" of  Ep  to  E$p. 

C.  The  Far  Scattered  Fields 

This  section  gives  the  derivation  of  the  formulas  for  the  far 
scattered  fields  of  a  ring  of  current  having  p,  $  and  z  components  of 
current  density.  The  far  fields  of  a  cylindrical  shell  of  current  are 
then  obtained  by  superposition  of  the  fields  of  the  elemental  rings  making 
up  the  shell.  1 

The  far  field  of  a  source  in  terms  of  cylindrical  source  components 
as  given  by  RichmondC28]  is  given  in  Eq.(34).  The  primed  coordinates 
are  the  source  coordinates  and  Fig.  8  defines  the  geometry  for  the  far 
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For  the  case  of  a  very  short,  thin  ring,  of  length  i,  thickness  t 
and  mean  radius  a,  the  integrations  over  p'  and  z'  can  be  replaced  by  a 
multiplication  by  (t -i).  In  addition,  if  we  replace  Jp ,  and  Jz  by 


field  calculation. 
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Fig.  8—Geometry  for  the  far  field  calculation. 


the  equivalent  current  of  Eqs.  (4),  we  find,  for  the  scattered  fields 
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where  k2  -  k  2  =  w2p  Je  -e)  . 
e  0  0 

The  field  expansion  of  Eqs.  (7)  and  (8)  is  used  in  Eq.  (35)  and, 
in  addition  we  let  cm  =  cos  (^r-z),  sm  =  sin(^-z)  and 
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in  order  to  find  (for  the  n-th  mode  in  $): 
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For  the  cases  considered  in  this  paper  the  incident  field  is  a 

7J 

plane  wave  at  <p  =  -  .  If  this  value  of  <j>  and  the  Fourier-Bessel  series 
2 

(Eq„  (19))  are  used  in  Eq.  (36),  the  final  result  for  the  far  scattered 
fields  of  a  ring  of  "equivalent  current  density"  can  be  shown  to  be: 
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The  scattering  cross  section  of  a  scattering  body  Is  defined  as 
follows: 


(38) 
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In  order  to  calculate  the  scattering  cross  section,  then  It  is 
simply  a  matter  of  calculating  the  far  scattered  field  and  then  applying 

t 

Eq.  (38).  It  should  be  noted  that  £in  Eq.  (37)  involves  a  factor  1/r 

p 

Which  eliminates  the  r  In  Eq.  (38).  E^ ,  of  course.  Is  the  field  in¬ 
cident  upon  the  body  and  in  this  case  is  taken  to  be  of  unit  magnitude 
so  that  |E^|2  =  1. 
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Fig.  9— Geometry  of  dielectric  ring  problem 
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(40) 

cont. 
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For  a  linearly  polarized  plane  wave  incident  along  the  ring  axis, 
only  the  n  =  1  mode  in  <P  is  needed.  This  situation  represents  the  ut¬ 
most  simplification  of  Eq.  (17).  Calculated  and  measured  backscatter- 
ing  data  are  compared  in  Fig.  10  for  several  polystyrene  rings  of  dif¬ 
ferent  radius  where  e  is  taken  to  be  2 . 54- j  0.00.  Each  ring  considered 

r 

had  a  geometrical  wall  cross  section  0.100"  x  0.100".  The  scattering 
cross  section  was  measured  twice  for  each  ring;  where  only  one  circle 
is  evident,  the  experimental  data  were  in  nearly  exact  agreement  for 
both  measurements.  The  agreement  between  the  calculated  and  measured 
results  is  seen  to  be  excellent. 

As  can  be  seen  from  Eq.  (16),  the  formulas  defining  the  coefficients 
in  the  system  of  linear  equations  are  independent  of  the  dielectric 
constant.  The  dielectric  constant  appears  only  in  the  solution  of  the 
set  of  simultaneous  equations  and  in  the  calculation  of  the  far  scattered 
fields;  these  two  operations  represent  relatively  efficient  computer 
operations.  (The  numerical  integrations  and  solution  of  the  simultaneous 
linear  equations  occupy  most  of  the  computation-time.  These  need  not 
be  repeated  for  each  new  permittivity  or  loss  tangent.)  Thus,  once  the 
basic  integrations  have  been  performed,  it  is  a  matter  requiring  rela¬ 
tively  little  additional  labor  to  calculate  the  effect  of  the  dielec¬ 
tric  constant  and  to  include  the  effect  of  the  loss  tangent  on  the  scat¬ 
tered  fields.  This  is  an  interesting  characteristic  of  the  integral 
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equation  approach  to  homogeneous  dielectric  scattering  problems. 
Fig.  11  shows  the  calculated  backscattering  cross  section  versus 
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Fig.  10— Backscattering  cross  section  versus  mean 
ring  radius  for  a  polystyrene  ring  with  a  plane, 
linearly  polarized  wave  Incident  on  axis. 
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dielectric  constant  for  a  plane  wave  incident  along  the  axis  of  a 
ring  with  an  outer  radius  of  0.5  inches.  The  effect  of  the  loss 
tangent  Is  also  included  parametrically. 

A  comparison  of  the  calculated  and  experimental  backs cattering 
cross  section  for  TE  or  TM  plane  waves  incident  at  an  arbitrary  in¬ 
cidence  angle  is  presented  next.  Four  polystyrene  rings  were  considered. 
Calculations  were  performed  for  several  incidence  angles,  out  to  nearly 
90°,  for  each  ring.  Both  TE  and  TM  backs cattering  patterns  were  meas¬ 
ured  for  all  four  rings.  Figures  12  through  19  show  the  comparison  of 
the  measured  and  calculated  data.  The  agreement  is  seen  to  be  excellent. 

A  portion  of  the  small  disagreement  in  the  TE  cases  may  arise 
from  scattering  from  the  support  strings;  the  strings  were  parallel  to 
the  incident  field  in  this  case.  For  the  TM  case,  the  incident  field 
was  normal  to  the  support  strings  and  the  string  scattering  was  not  so 
noticeable. 

The  number  of  modes  in  <p  required  for  a  particular  ring  scattering 
problem  depend  on  the  ring  radius  and  the  propagation  angle  of  the  in¬ 
cident  field,  i.e.,  the  nunfcer  of  incident  modes  which  have  significant 
magnitude.  In  the  case  discussed  here,  both  the  incident  modes  and  the 
far-fleld  scattered  modes  contain  factors  Jn+^(ua)  or  Jn_^(ua);  the  prob¬ 
lem  of  determining  the  necessary  number  of  modes  reduces  to  an  examin¬ 
ation  of  the  Bessel  functions  of  the  first  kind.  In  the  examples  pre¬ 
sented  here,  the  incident  field  modes  became  negligibly  small  for  n  _>  10  ; 
therefore  only  ten  $-modes  were  required.  For  more  complicated  incident 
fields  and  for  larger  rings,  more  modes  may  be  required;  the  exact  number 
would  of  course  depend  on  the  accuracy  desired. 


POLYSTYRENE  RINO 
<  -  2.54  -  jO.OO 
FREQUENCY  -  6.030Hz 


44 


Fig.  14— Plot  of  backscattering  cross  section  versus  incidence 
angle  for  a  TM  plane  wave  incident  on  a  dielectric  ring. 
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Fig.  16--Plot  of  backs cattering  cross  section  versus  incidence 
angle  for  a  TM  plane  wave  incident  on  a  dielectric  ring. 
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Fig.  18--P!ot  of  backscattering  cross  section  versus  incidence 
angle  for  a  TM  plane  wave  incident  on  a  dielectric  ring. 
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Fig.  19— Plot  of  backscattering  cross  section  versus  incidence 
angle  for  a  TE  plane  wave  incident  on  a  dielectric  ring. 
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It  should  be  noted  that  the  only  Inherent  assumption  in  the  ring 
scattering  calculation  is  that  the  thickness  and  length  are  sufficiently 
small  compared  with  a  wavelength  that  field  variations  with  p  and  2  can 
be  neglected.  Aside  from  practical  problems  such  as  lengthy  Integration 
times  and  summations  over  many  $-modes,  there  is  nothing  preventing  ap¬ 
plication  of  this  technique  to  rings  with  larger  radii  If  the  above  re¬ 
strictions  are  kept  in  mind. 

B.  Dielectric  Tube  Scattering 

The  final  situation  considered  is  that  of  a  homogeneous  thln-wall 
cylindrical  shell  with  a  linearly  polarized  plane  wave  Incident  along 
the  shell  axis.  Shells  up  to  one  wavelength  long  are  considered  and 
several  terms  of  the  z-expanslon  are  required  for  a  good  solution  for 
the  field  In  the  tube.  Choice  of  on-axis  Incidence  requires  only  the 
n  ■  1  mode  In  4  and  the  amount  of  computer  storage  required  Is  consid¬ 
erably  reduced  from  the  arbitrary  Incidence  calculation.  The  geometry 
of  the  problem  is  shown  in  Fig.  20. 

In  this  calculation,  all  the  coefficients  given  in  Eq.  (16)  were 
retained  In  the  system  of  equations.  Although  rapid  access  computer 
storage  restrictions  limited  the  maximum  z-mode  Index  to  two,  this  was 
found  to  be  adequate  for  the  cases  considered.  For  tubes  much  larger 
than  one  wavelength,  however,  more  modes  would  be  needed. 

The  calculations  and  experiments  were  performed  for  polystyrene 
tubes  with  an  outer  radius  of  0.500"  and  a  wall  thickness  of  0.100". 

The  frequency  was  6.23  GHz,  the  dielectric  constant  used  in  the  calcu- 
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Fig.  20--Plane  wave  incident  along  dielectric  tube  axis. 


lations  was  er  =  2.54  -  jO.OO  and  the  tube  lengths  ranged  from  0.100" 
to  1.960". 

The  experimental  data,  the  averaged  data  and  the  calculated  back 
scattering  cross  section  are  given  in  Table  I.  A  comparison  of  the 
caculated  and  measured  scattering  cross  section  is  shown  in  Fig.  21. 

The  agreement  between  the  calculated  and  measured  cross  sections 
again  is  excellent.  The  increasing  disagreement  for  the  larger  tubes 
undoubtedly  arises  from  the  three-mode  limitation  in  the  z-expansion. 
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CONCLUSIONS 

Accurate  calculation  of  radome  boresiqht  error  is  a  challenging 
task.  Rigorous  series  solutions  of  the  boundary-value  problem  are 
usually  not  practical  for  the  general  radome  shell  and  the  accuracy 

I 

of  present-day  approximate  methods  of  radome  analysis  is  generally 
not  known.  An  essentially  rigorous  technique  for  application  to 
radome  scattering  is  the  integral  equation  method.  The  first  step  in 
the  development  of  an  integral  equation  analysis  of  radome  shells  is 
discussed  in  this  report,  namely  the  calculation  of  the  backscattering 
from  thin-wall,  finite-length,  circular  cylindrical  shells. 

The  unknown  function  in  the  integral  equation  is  the  field  within 
the  dielectric  shell.  The  integral  equation  derivation  and  the  method 
of  solution  of  the  equation  by  the  point  matching  technique  are  dis¬ 
cussed.  The  field  within  the  dielectric  shell  is  found  In  terms  of 
expansion  coefficients  in  a  double  Fourier  series  over  the  shell.  The 
far  scattered  fields  are  determined  by  using  the  calculated  total  field 
and  the  equivalent  current  concept. 

The  technique  is  applied  to  calculation  of  the  backscattering  of 
a  linearly  polarized  Incident  wave  by  dielectric  rings  and  tubes.  A 
comparison  of  the  calculated  and  measured  backscattering  was  made  for 
the  following  three  particular  cases: 

(1)  a  thin,  short,  homogeneous,  dielectric  ring  for  axial 
Incidence 
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(2)  a  thin,  short,  homogeneous,  dielectric  ring  for  arbitrary 
incidence  with  TE  and  TM  polarization 

(3)  a  thin-wall  cylindrical  shell  for  axial  incidence. 

The  dielectric  material  for  each  of  the  above  cases  is  polystyrene. 
Agreement  between  the  calculated  and  measured  results  is  excellent  in 
all  three  cases. 

The  general  numerical  integration  procedures,  the  far  field  in¬ 
tegration  formulas,  the  experimental  methods  and  the  computer  programs 
used  are  all  discussed. 

In  addition,  a  useful  expression  is  derived  for  the  electric  field 
at  the  center  of  a  rectangular  parallelepiped  of  current  density.  This 
expression  is  needed  when  integrating  throughout  a  volume  of  source 
current  density  to  calculate  the  field  at  a  point  within  the  source 
region. 

The  numerical  technique  used  seems  to  be  particularly  suitable  for 
treatment  of  arbitrary  shaped  shells  of  revolution  which  could  be  con¬ 
sidered  to  be  made  up  of  a  large  number  of  dielectric  rings.  Longitu¬ 
dinally  inhomogeneous  bodies  are  also  amenable  to  analysis  by  the  ring- 
subdivision  method;  each  ring  would  be  homogeneous,  of  course,  so  that 
this  approach  would  give  a  piecewise  uniform  approximation  to  the  in- 
horoogeneity  being  considered. 

In  summary,  excellent  results  have  been  obtained  by  applying  an  in¬ 
tegral  equation  technique  to  calculation  of  the  plane  wave  scattering  from 
thin-wall  dielectric  cylindrical  shells.  This  work  Is  significant  to  the 
general  radome  scattering  problem  for  two  reasons.  First,  the  excellent 
theoretical  and  experimental  agreement  implies  a  high  degree  of  accuracy 
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obtainable  with  the  nethod.  Second.  .  portable  node  of  eviction  of 
the  Integral  e, nation  technlgue  has  been  established  and  this  approach 
can  be  extended  to  »ore  co^e*  d1e,ectr1c  shell  configurations  and  to 

more  general  fields. 


APPENDIX  A 


EXPERIMENTAL  METHODS  AND  EQUIPMENT 
Measurement  of  the  backscattering  cross  section  of  the  rings  and 
tubes  was  accomplished  using  a  one  horn  monostatic  backscattering 
cross  ection  system  operating  in  a  large  microwave  darkroom.  A  block 
diagram  of  the  measurement  system  is  shown  in  Fig.  22. 

A  measurement  is  made  by  following  the  steps  outlined.  The  three- 
stub  tuner  is  adjusted  to  create  a  load  mismatch  precisely  balancing  out 
the  background  return  from  the  darkroom  when  no  target  is  present.  A 
reference  cross  section  level  is  then  established  by  using  a  standard 
sphere  with  known  cross  section  as  the  target  and  recording  the  back- 
scattering  return.  The  sphere  is  then  removed,  the  target  is  placed 
in  position  and  the  target  return  is  recorded.  The  absolute  value  of 
the  target  cross  section  is  then  found  by  measuring  the  relative  levels 
of  the  target  and  sphere  and  calculating  the  target  cross  section 
directly  from  the  known  sphere  cross  section.  In  all  cases  the  reference 
sphere  is  chosen  to  provide  a  reference  level  of  the  same  order  of  mag¬ 
nitude  as  the  target  cross  section. 

The  linearity  of  the  system  was  checked  by  comparing  the  measured 
levels  of  two  reference  spheres.  For  the  range  of  cross  sections  con¬ 
sidered,  the  relative  sphere  levels  were  within  -0 .1  dB  of  their  known 
relative  values.  This  is  considered  to  be  very  good.  The  wide  range 
linearity  (30  dB  range)  was  also  checked  using  a  precision  attenuator; 
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Ftg.  22— Block  diagram  of  moncst.atic  cross  section 
measurement  app>r-t.us. 
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this  check  gave  relative  levels  within  O.ldB  of  the  attenuator  values. 

For  the  case  of  the  on-axis  scattering  of  the  thinnest,  shortest 
rings  (Fig.  10),  the  cross  section  levels  were  somewhat  low  to  obtain 
good  "far  zone"  scattering  results.  As  pointed  out  by  Kouyoumjian  and 
Peters, T29]  however,  accurate  cross  section  measurements  do  not  require 
that  the  target  be  in  the  far  zone  of  the  antenna  but  only  that  the 
Incident  fields  have  plane  wave  properties  in  the  vicinity  of  the  tar¬ 
get  and  that  the  antenna-target  interaction  be  very  small. 

The  thin,  short  rings  were  measured  by  suspending  them  on  strings 
in  the  horn  mouth  with  the  center  of  the  rings  on  the  horn  axis. 
KouyoumjianC30]  shows  that  in  the  neighborhood  of  maximum  points  along 
a  horn  axis,  the  field  is  nearly  a  uniform  plane  wave.  The  axial  field 
variation  of  the  measurement  horn  was  probed  using  a  0.25"  sphere.  For 
the  thin  rings,  the  axial  incident  field  variation  over  the  rings  was 
found  to  be  less  than  0.5  dB  and  the  horn-mouth  location  was  found  to  be 
very  close  to  a  maximum  point.  For  small  rings  and  small  cross  sections, 
antenna-ring  interaction  can  undoubtedly  be  neglected.  Thus,  reasonably 
good  plane  wave  scattering  measurements  are  expected.  No  quantitative 
estimate  of  the  errors  has  been  made,  however. 

Measurement  of  the  ring  patterns  and  the  tube  cross  sections  Im¬ 
poses  a  larger  axial  extent  of  the  target  than  can  be  tolerated  for  a 
horn-mouth  measurement.  The  larger  cross  section,  however,  allowed 
use  of  the  more  standard  far-zone  cross  section  measurement.  Observed 
repeatability  and  linearity  were  again  very  good. 


APPENDIX  B 
COMPUTER  PROGRAMS 


A  discussion  of  the  computer  programs  used  to  calculate  the  ring 
and  tube  fields  and  backscattering  cross  section  is  presented  in  this 
appendix.  The  programs  are  written  in  the  Fortran  IV  language  for 
processing  on  the  Ohio  State  University  IBM  7094  computer. 

Flow  charts  outlining  the  computations  are  given  for  all  the  pro¬ 
grams  used.  Definition  of  the  input  and  calling  parameters  and  a  com¬ 
plete  statement  listing  of  the  programs  are  also  included. 

A.  Arbitrary  Plane  Wave  Incidence  on  a  Dielectric  Ring 

The  program  flow  chart  is  shown  in  Fig.  23  and  the  statement  list 
ing  of  the  computer  program  is  given  in  Fig.  24.  The  input  data  card 

variables  and  the  subroutine  calling  paramenters  are  defined  below. 

(1)  Cards  needed  once  for  each  run: 


Card  -Number  1  Format  (110,  7F  10.5) 

Columns  Quantity  Description 

j.10  NCASE  NCASE  *  the  number  of 

different  rings  to  be 
calculated. 


(2)  Cards  needed  once  for  each  different  ring: 

Card  Nunter  2  Format  (110, 7F  10.5) 

Columns  Quantity  Description 


1-10 

(right-adjusted)  NX 

11-20  AI 


NX  «  the  maximum  ^-mode 
index  to  be  used 

AI  =  the  mean  ring 
radius  in  inches. 


d 
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1 

Columns 

Quantity 

Description 

1 

i 

1 

i 

21-30 

DI 

DI  =  the  radial  thickness 
of  the  ring  in  inches. 

31-40 

WI 

WI  =  the  width  of  the 
ring  in  the  z-di recti  on 
in  inches. 

41-50 

FREQ 

FREQ  =  the  frequency  in 

1 

. 

gigahertz. 

I 

51-60 

Re(ER) 

ER  *  the  complex  relative 
dielectric  constant  of 

4 

61-70 

Im(ER) 

the  ring  material. 

] 

Card  Number  3 

Format  (I10.7F  10.5) 

• 

1 

Columns 

1-10 

Quanity 

Description 

(right-adjusted) 

NANG 

NANG  =  the  number  of  in¬ 
cidence  angles  to  be 

! 

• 

calculated  for  the  ring. 

f 

11-20 

REF 

REF  =  the  reference 
cross  section  for  calcu¬ 
lation  of  the  theoreti¬ 
cal  cross  section  in  dB. 

• 

(3)  Cards  needed 

once  for  each  different  incidence  angle  for  each  ring. 

» 

Card  Number  4 

Format  (8F  10.6) 

* 

Col  umns 

Quantity 

Description 

1-10 

11-20 

21-30 

B(I),  1=1,10 

B(I)  *  the  I-th  order 
Bessel  function  of  the 
first  kind.  Used  in 

etc. 

calculation  of  the  in¬ 
cident  field  Fourier  co¬ 
efficients  for  a  partic¬ 
ular  incidence  anqle. 

i 

Card  Number  5 

Format  (8F  10.6) 

? 

i 

Columns 

Quantity 

Description 

• 

1-10 

TH 

TH  *  the  incidence  angle 
for  one  calculation  for 

one  ring. 
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The  calling  parameters  for  the  NRCELL  subroutine  are  defined  as 
fol lows: 

M  =  the  number  of  subcells  along  each  side  of  the  larger  cell 
(a  total  of  subcells) 

RC,PC,ZC  are  the  0,4,2  coordinates  respectively  of  the  larger 
(near)  cell 
AL  =  O.d 

NX1  *  the  maximum  4-mode  index  to  be  used  plus  1. 

XM,YM,ZM  =  the  rectangular  coordinates  of  the  matching  point  for 

which  the  NRCELL  subroutine  was  called. 

« 

AA,BB,CC  *  the  incremental  dimensions  of  the  large  (near)  cell  in 
the  p,4  and  z  directions  respectively. 


Main  Program 


Add  the  singular  cell  contribution  to 
the  appropriate  coefficients.  Do  for 
all  4-modes. 


Punch  out  (on  cards)  the  coefficients 
in  the  system  of  equations.  Do  for  all 
4-modes . 


400  continue 


Do  800  N  =  1,  NANG 


Read  Bessel  functions  and  incidence 
angle. 


1002  continue 


Calculate  incident  field  Fourier 
coefficients . 


35  continue 


Solve  system  of  equations  for  Fourier 
coefficients  of  the  total  field. 


500  continue 


Calc,  the  far  scattered  fields  for  all 
i-modes . 


600  continue 


Calc,  the  total  scattered  field  as  the 
sum  of  the  fields  due  to  each  of  the 
•-modes . 


Calculate  and  write  out  the  backscattering 
cross-section 


800  continue 


NRCELl  SUBROUTINE 


Calling  parameters  (M,RC,PC,ZC,  AL,  NX1, 
XM,  YM,ZM,AA,BB,CC) 


Initialize  integrals  to  zero 


50  continue 


Divide  the  cell  into  subcells  and  cal¬ 
culate  the  total  cell  integral  as  the  sum 
of  the  contributions  of  each  of  the 
subcells . 


r>  r»  r»  r»  on  fvnftftftn 
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VANOOEREN.  R.  J08N  FGJ921  3 

•EXECUTE  IB JOB 

•IB JOB 

•IBFTC  MAIN  NOOECK 

1  FORMAT  Cll0.7Fl0.5) 

2  FORMAT  ( I5X.2IS) 

22  FORMAT  (5x.aFl4.BI 

3  FORMAT  I  4X  .  21  H  REST  OF  RING  RESULT S/4 X  *  1 4F8.4 ) 

4  FORMAT 16X.21H  SINGULAR  CELL  RESULT/1  OX . 4F8 .4 1 
7  FORMAT  (5F13.8) 

15  FORMAT  (BF10.6I 

COMPLEX  A | ( 10 ) . A3 ( 1 0  > . A4< I  0 ) . A6< 10 ) . A8  t |0 ) 

COMPLEX  AT< 10) .BT1 |0 ) . API tO I .8P< 10). ALl IOI.BLI 10) 

COMPLEX  AlT<lO).eiT(lO).A[P(lO).BtP(10>.AILflO>.BILUO> 

COMPLEX  CJ«Ct .ETH1 10)i EPS < I  0 ) .EThT.EPST 
COMPLEX  XP.P.O 
COMPLEX  det.er 

COMPLEX  VI  I  1 0 ) «V2< 1 0) .V3< I  0  >  >V4 1 1 6) . V5< 10) .V6< 10) »V81 10 ) 
COMPLEX  V9<10).VI1 ll0).V16(10> 

complex  pi < i O) .P2i i o  >  «P3t i o  >  «pa  t io> .psi 10 ) .P6c iO) «p71 10 ) »pbi io) 
COMPLEX  P9<10).PlO<lO>.Ptl(lO).Pl2110>.PI3(IO>,P14t|0).P15<IO> 
COMPLEX  P 1 6 ( 1 0 ) 

COMMON  PI .P2.P3.P4.P5.P6.P7.P8.P9.Pl0.Pl 1 .P12.P13.Pl4.PI5.Pl6 
COMPLEX  AITMI 101.8ITMI 10 ) . A1PMI 10 ).BIPM« 10). AlLMI 10).BILM( 10 ) 
COMPLEX  ATM ( 1 0 ) .BTM 110). APM  HO) *BPM 1 1 0 ) . ALM 110) .BLM 110) 

COMPLEX  ETHMi 10).EPSM1 10) .ETHMT.EPSMT 
DIMENSION  B< 10) 

COMPLEX  EE 
COMPLEX  XPP.XPO 
RE AO  <5.1 )  NCASE 

C  NCASE  IS  THE  NUMBER  OF  RING  CASES  TO  BE  RUN. 

00  BOt  NCAS«|.NCASE 
RE AO  15.1)  NX  « A | «DI .FI »FREQ»  ER 
NX  IS  THE  MAXIMUM  PH I -MODE  INDEX. 

At*  MEAN  RADIUS  OF  THE  RING  IN  INCHES. 

Dl»  THICKNESS  OF  RING  IN  RADIAL  DIRECTION  IINCEHES). 

WU«WIDTH  OF  RING  IN  Z-DIRECTION  1  INCHES). 

ER  a  RELATIVE  OIELECTRIC  CONSTANT  OF  RING. 

FREO  a  FREQUENCY  IN  GHZ. 

READ  15*1)  NANG. REF 

NANG  IS  THE  NUMBER  OF  INCIDENCE  ANGLES  FOR  A  GIVEN  RING. 

REF  IS  THE  REFERENCE  CROSS-SECTION  FOR  THE  PARTICULAR  RING. 
AU«0.0 
nxi«nx+i 

Y00«  U.8028/FPEO 
Vila  Y00/SQHT1ER) 

AaAI/YOO 
D«0I/YO0 
WaWI/YOO 
WAO*6.283l9/VOO 
WA»6.28 )I9/Y1  1 

A  IS  THE  MEAN  RADIUS  OF  THE  RING. 

D  IS  THE  THICKNESS  IN  THE  RADIAL  DIRECTION  OF  THE  RING. 

W  IS  THE  WIDTH  OF  THE  RING  IN  THE  Z-DIRECTION. 

ALL  THE  ABOVE  ARE  IN  TERMS  OF  FREE  SPACE  WAVELENGTHS 
00  32  J«1.10 
VI «J)f<0..0. 1 
V2I J)»<0.,0. ) 

V3(J)f(0..0. ) 

V41 J>«(0..0. ) 

VS  1 J)»  <0. . 3. 1 
V6<J)«<0..0. I 
V8< J)«10. ,0. ) 

V91 J)a(O.«0. ) 

Fig.  24--Ring  scattering  computer  program 
statement  listing. 


VII  <  J1K0..0.  » 

V16(J1«(0..0. I 
ATI J)*<0..0. ) 

8T(J1*<0..0. ) 

AP(J1*< 0..0. > 

8P( JI*<0..0. > 

AL<J1I:<0..0« ) 

8L(J1=(0..0. ) 

AIT( J1*(0.,0. » 

BITIJIXOmO.  ) 

AIP< J1«<0..0. ) 

8IP< Jl*<0..0. ) 

AIL< J1*<0. <0. > 

8JL< J1*C0..0. ) 

ETH< J>*<0..0. ) 

EPS(U1*<0..0. 1 
ETHMC Jl*<0..0.1 

epsm<ji*<o..o,i 
32  CONTINUE 

ETHT«(0.«0.'l 
EPST*(0..0.1 
ETHMT*<0.,0. ) 

EPSMTXO.'O.  ) 

00  403  I NO* I • NXI 
READ  (5.16)  J.AI.FREO 

READ  (3.171  A| ( J) «A3( J) <A4( JI.A6IJ). A8( J> 

403  CONTINUE 
00  TO  401 

1000  CONTINUE 

N*  ( 3. 1 4159* A  1/0.0 1 
OGL*0.01/A 

C  THE  CELL  FROM  -DEL  TO  TOEU  IS  ACCOUNTED  FOR  ANALYTICALLY. 

IF  (N.LT.30)  N«30 
T«*N 

OPHI* <3.1  4139-DEL 1/T6 
OV*W*D*A*DPHI/6.28319 
WTT * ( (W/2, 1**31*0.666667 

1001  CONTINUE 

OO  200  K*20.N 
T 1  mK 

PH»OEL +0PHI/’2.  +  TI  *0  PH  I 
-  C*COS(PHl 
S*SI NCPH) 

X»A*( 1 .-Cl 

Y*-A*S 

Z«0.0 

R*SORT ( X*X*Y* Y+Z*Z l 

AR*6.2BJI9*R 

ARS»AR**2 

R3«R**3 

R9«R**S 

PR*  (-1  .+ARS1/R3 

PI •— AR/R3 

P*CMPLX(PR.P|  1 

OR* (3.— ARS 1 /R5 

Q I  *3.*AR/'R3 

0*CMPLX(0R.0I I 

XP*CMPLX(COS< ARI.-SIN(AR) 1 

XPP*XP*P 

XP0*XP*0 

00  60  NN*1 .NXI 

T1*NN-1 

4N«SIN(T1*PH1 

CN*COSt  T1 *PH 1 

Fig.  24—  (continued). 
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CNO»CN*OV 
SNO«SN*0V 

VI  (NN)«VI  (NN)+C»CNO«XPP  ' 

V2 (NN )  «V 2  (NN  )  +S*SN0*XPP 
V3(NN).V3<NN)+X«X*C*CN0*xPQ 
V4(NN)«V4  (NN).X*X*S*SN0«XPQ 
V5 (NN ) »V5  <NN ) .X*Y»S*CN0*XP0 
V6  (NN  )  *V6  (NNM-X*V*C*SNO*XPQ 
V8(NN)*V8<NN).Y*Y*S*SN9»XPQ 
V9(NN)*V9(NNJ+Y*Y#C*CND*XPO 
VII (NN)iVtl (NNI*CNO*xPP 
V 1 6 ( NN ) *V 1 6 ( NN ) ♦WTT  *CNO*XPQ 
60  CONTINUE 

200  CONTINUE 
00  201  J«1 .NX! 

WRITE  (6.1)  J 

WRITE  (6« 3 )  VI ( Jl 4 V2 I J) . V3 1 J > «V4 I J > . V5C J ) . V6< J ) «V8 ( J) . V9( J) . VI I ( J) 

201  CONTINUE 
XM«A 
YH«0 . 0  ) 

ZM*0.0  1 
4*«0 

BB«S.*0PH1 
CC«W 

OOWN  TO  250  IS  THE  CALCULATION  FOR  THE  4  CELLS  NEAREST  THE  *** 
***  SINGULAR  POINT  (  THE  OBSERVATION  POINT). 

DO  250  J«1 .4 
TI«J-I 
OPH*5»*OPHl 
PH I «DEL*0PH/2 . +T 1 *OPH 
C-COSIPH! ) 

S*SIN(PHI I 
XC»A*C 
YC«A*S 
2C«0.0 

IF  (J.EO.l)  M«20 

IF  (J.E0.2I  H«lO 

IF  ( J.EQ.3 )  M»8 

IF  (J.E0.4)  M»6 

CALL  NR '"ELL  (  M.  A  .PH  I  .  ZC  .  AL  .  NX  1  ,  XM.  VM  .  ZM  .  A  A  ,BB.  CC  ) 

00  202  JJ«I.NX! 

WRITE  (6.2)  J.J J 

WRITE  (6.22)  Pl ( JJ) .P2( JJ) «P3( JJI.P4 ( JJ ) «P5 ( JJ ) «P6 ( JJ ) .PT ( JJ ) .PB ( J 
2J) .P9( JJ) .PlO( JJ) .P| 1 ( JJ ) »P| 2  I JJ  I  «P| 3( JJ  >  «P) 4 ( JJ ) «Pl 5( JJ ) .Pi 6( JJ ) 
VI < JJ)«V1 ( JJ)*P1 ( JJ) 

V2 ( JJ ) *V2 ( J J ) +P2 ( J J ) 

V3 ( JJ  > ■ V3 ( J J ) 4P3 ( J J ) 

V4  (  J  J )  *  V4  ( J  J  )  4P4  (  J  J  ) 

VS( JJ ) »V5 ( JJ ) +P5 ( J J ) 

V6 ( JJ ) « V6 ( JJ ) *P6 ( J J ) 

VB ( J J ) «V8 ( JJ ) .PS ( J J ) 

V9 ( J J ) » V9  <  J J ) APR ( J J ) 

V 1 1 ( JJ ) *V 1 1 ( J J)*P II ( JJ) 

V 1 6  ( JJ )  «V  1 6  (  J  J  )4P1 6  ( J  J) 

WRITE  (6.2)  JJ 

WRITE  (6.22)  VI (JJ)«V2(JJ>.V3(JJ>*V4(JJ).VS(JJ).V6(JJ)»VS{JJ>.V9(J 
2  J  )  .  VI  t  ( JJ  1 .  V 1 6  ( J  J ) 

202  CONTINUE 
260  CONTINUE 

c  FROM  HERE  TO  300  IS  CALCULATED  THE  SINGULAR  CELL  CONTRIBUTION. 

AI-0/2. 

B I »A»OEL 
C!«W/2. 

SQ. SORT (A  1 *A | *81 *BI *C1 *C I  ) 


Fig.  24--{continued) 
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T2*2./3.l4159 

4T«T2*/ATAN<  <AJ»CI  }/<Bl*SOl H-ATANt  <A1*8I )/(Cl*SO) ) 1-1 « 
SP*T2*(ATAN< (B!*A1 )/(Cl*SQ) l+ATANI (B1*CI i/(Al*S01 1 1-1. 
5L*T2#(ATAN< (Cl *81 !✓< A|*S01 l+ATANI <C1*AI }/CBt*S0! ) >-t . 

300  CONTINUE 

c  ER  IS  THE  RELATIVE  DIELECTRIC  CONSTANT  OF  THE  DIELECTRIC  RING. 

00  400  J*t .NX1 
At  <J!»VI  < J)+V3(UM-V5I  JlfST 
A3( J)«-V2( Jl-V4< JI+V6I J) 

A4< J1*V2< Jl+V6< JI+VS( Jl 

A6( J1*VI < JI*V9< JI-V5I J1+5P  . . . 

A8( J)»V1 1 < J)+V16( JI+SL 

WRITE  (6.S)  J.A1 ( U 1 «A3(J1 « A4 (J)*A6(J).A8(J) 

16  FORMAT  (H0.7F10.6) 

17  FORMAT  (8F10.6) 

PUNCH  16.  J. At .FREQ 

PUNCH  17.  A 1 ( J 1 . A3 ( J I . A4  <  J I .  A6 ( U  >  * AB ( J 1 

400  CONTINUE 

401  CONTINUE 

DO  BOO  N1 1*1 .NANG 

READ  <5.151  !(8<J).J*1.101 

C  THE  B.S  ARE  BESSEL  FUNCTIONS  USED  IN  THE  INCIDENT  FIELD  CALC. 

_READ  (5.151  TH  . 

1 002  CONTINUE 

ST*SIN(J. 14159- TH! 

CT«COS ( J • 1 4 1 59— TH 1 
TM0»TH*57.3 

C  DOWN  TO  35  IS  CALC.  OF  THE  INCIDENT  FIELD  COEFFICIENTS. 

_  A 1 T  (  1  1X0. .0.  1 

B IT ( 1 1* (O. .0.1 
A!P(1  1«<-MPLX(0..-B(21  1 
BIP(1 1*10. .0. 1 
A ITM  < 1  !*-CT*A IP ( 1 1 
BITM< 1 1X0. .0.} 

AlPMd  IXO.iO.) 

BtPMd  1X0. .0.1 

AILMd  1*-ST*CMPLX(B(1  1.0.1 

BtLMd  !«<0..0.  | 

DO  35  JX .4 
J1»2*J 
JlPl-Jl+1 
J1M1-J1-1 

-AIT  (  J 1  1-CMPLXT  <S< JIM1  1+B(  JIPI  11.0.  1 
BIT  (  J1  IX0..0.  I 
A|P( J1 1»<0..0.1 

BIP( J1 1 *CMPLX ( <  — B ( Ul Ml !7B( JlP| 1 1.0. I 
AITMIJ1 1*<0..0.1 

.  B I TM ( J 1 l«-CT*BIP< J1 1 

A I PM ( J I  !«CT*A|T(JI I 
B1PM(Jt  XIOxO.I 
AILMIJl  l«<0. .0.1 

B ILM ( J 1 |.-CMPLX(0..2.*ST*B(J1 1 1 
J2«2*J+1 
“  '  J2Pl*J241 

J2M1«J2-I 
A 1 T  (  J2 1  *  <  0  .  .  O  .  ) 

BIT ( J2 1 *CMPLX (0. . (B ( J2M1  |♦B(U^PI I  1  I 

A1P( J2I«CMPLX(0..(B(J2M1 ]-B(J2P| 111 

B1P(  J21X0..0.1 

A|TM( J2 I*— CT*A|P<  J2 1 

BITMIJ21X0..0.1 

A|PM(J2  1X0.  .0.  1 

BIPM(J2!«CT*BIT(J21 

AILMI J2!*-CMPlX(2.*ST»B( J21 .0.1 


Fig.  24— (continued) 
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8 ! LM  t  J2  X  <  0  .  « 0  .  ) 

35  CONTINUE 

FOLLOWING  IS  The  solutine  fob  the  foubier  coefp.  op  the  fields. 
EE»ER-1 . 

DO  500  Jaj ,NX1  _  ... 

OETalj .-EE*Al <  J » >  *  < l.-EE*A6( J) XEE*A3<  J  XEE*A4  <  J  > 
AT<JX<AITIJ)*I1.-EE*A6(J))+BIP(J )*EF**3  <  J >  >/DET 
8P< JX  <BIP< J)*< I .-EE*A1 < J) )-AtT< J)*EE*A4< J) l/DET 
API  JX<A1P<J)*<1.-£E*AI  <J>)+8IT(J)*EE*A4<J)  >/OET 
BT  <JX<OIT<JX<1 .-EE*A6< J) XA IP ( J >*EE* A3 < J > ) /DET 

AL  IJXIO..O.  )  ......  ...  _ 

8L<  JXIO..O.  > 

ATM < J 1 • < A  I TM < J »* < 1 . -EE* A6 < J » > +B I PM  I J ) *EE*A3 < J > | ✓DET 
BPM<JX<8IPM<  J)*<  |  .-EE*A1  <  J)  XA I TM I J > *EE*A4  < J )  >/DET 
APM<JX<A1PM<  J)*< I.-EE*A1 <J) > .8 I TM I J >*EE* A4 < J ) ) /DET 
BTM<JX<8ITM<  J>*<  I  .-EE*A6< J)  )-A  I  PM  <  J  )*EE*A3  <  J)  >/DET 

ALM  <  J  X  A  IlM  <J)/<1«-EE*A8<J)  I  . . .  . 

BLM  <  JXB1LM  <  J  >  /  <  1  .-EE* A Q  <  J  )  ) 

8  FORMAT  <SX. I  t0./.b:<*6F15.8./«5X.6Fl5.8»/) 

WRITE  <6. 18) 

WRITE  <6.61  J.AIT<J).SIT<J).AIP<J),BIP<J).A1L<J>.BIL<J> 

WRITE  <6.0)  J.AT<JxaT<J).AP<Jl.BP<J».AL<  JXBL<J> 

WRITE  <6. 19)  _ 

WRITE  <6.8  )  J.A1TMI  J  )  .B I  TM  <  J  )  ,  A  I  PM  I J  )  ,  B  I  PM  <  J  >  .  A  lLM  <  J )  ,B  Il.M<  J  I 
WRITE <6. 8)  J.ATM1 J).8TM<J).APM< Jl.BPMI J).ALM< J).BLM(J) 

500  CONTINUE 

THE  FOLLOWING  IS  CALC.  OF  THE  SCATTERED  FIELDS. 

CJXO.il.) 

ETH  < I )a-2.*B<2)*CT*CJ*AT<t 1+2. *8(1 )*ST#AL<I >  _ 

EPS<Vx-z.*a<2)*cj*AP<n 

ETHMIl X-2.*S<2)*CT*CJ*ATM< J X2.*B<1 )*ST*ALM<J ) 

EPSM<1 )«-2.*B<2)*CJ*APM< i ) 

OO  6 00  J*| o 
JI »2* J 

IF  <JI  .EQ.2.0R.J1.EQ.6.0R.J1  .EQ.IO)  Cl»-CJ  . 

IF  < J1.E0.4.OR. Jl.e0.6l  CXCJ 

JM»JI-1 

JPaJXI 

JMMaJM-t 

jPPaJP.l  • 

5I<M>)«J 

CTH< JP)aSI*CI*< <  <-AT<  JPXBP<  JP) 1*8 < JPP  > ♦ < AT  I JP XBP < JP) ) *B < JI  ) )*CT. 
22 • *AL I JP ) *  < “C  J)  *B  <  JP  >  *ST ) 

EPS  <  JP  XS 1  *C|  *  <  <BT<  JP)-A»  <  JP  )  )*B  <  JPP  )  +  <  8T  <  JP  XAP<  JP  )  »*B<  JI  )  ) 

WRITE  <6. 9) 

WRITE  <6.8 )  JP.ETH< JP) .EPS < JP) 

9  FORMAT  <5X.I8HJP.ETH(JP),EPS<JPI1  .  _ _ 

"  '  gTHMI JP)«SI*CI*<  <  <-ATM<JP)-BPM<JP ) )*8< JPP )+<ATM<JP)-BPM<JP> )*B< JI ) 

2>*CT*2.*ALM< JP)«<-CJ)*8< JP>*ST) 

EPSM  <  JP  )  ■  S  1  *C  1  *  <  <BTM<JPXAPM<  JP  )  >*B  <  JPP  X  (BTM  <  JP  XAPM<  JP  )  )»8  <  JI  >  ) 
WRITE  <6. 900) 

900  FORMAT  <5X.20HJP.ETHM< JP).EPSM< JP) ) 

WRITE  <6. 8)  JP,ETHM<JP),EPSM<JP)  .  _  . 

J8*  <2*J  X I 

JMaJ2-| 

jPa J2.1 
JMMaJM-t 
JPP. UP. 1 

IF  < J2.E0. I .0B.J2.EQ.5.0R. J2.E0.9)  C2»-1.0 
IF  < J2.E0.3.0B. J2.EQ.7)  C2«l«0 

J2Pa< J2.< \/Z 

SX<-l.)**J2P 

ETH< JP ) aSI*C2*<  <  <8T  <  JPXAPl  JP)  1 *M<  JPP )♦ < -6T  <  JP )-AP <  JP )> *8  I J2  >  >*CT ♦ 
22.*BL  <JP)*CJ*8< JP)*ST) 


Fig.  24--(continued) 
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EPS(  JP)»Sl*C2*(  (AT ( JP)+8P( JP)  )*B(JPP  >  +  I  AT  I  JP)-8P( JP) )4B( J2>  ) 

ETHM( JP)=S I *C2*< ( <8TM( JP)-APM( JP  1 )*B( JPP )  + (-BTMI JP l-APM  t JB ) ) *6t J2  » 
?)  4CT+2.*0L'->(  JP)*CJ*9(  JP)*ST  ) 

EPSMI JPl=S(*C2*( ( ATy ( JP)+BPM( JP) ) Ag < JPP  > +  < ATM ( JP )-BPM ( JP ) )4B( J2> ) 
WRITE  (6.10) 

10  FORMAT  (5X.ieHjP.ETH(jP>.EPS(JP) ) 

WRITE  (6.8)  JP,ETH( JP).EPS(JP) 

WRITE  (6.900) 

WRITE  (6.8)  JP.ETHM(JP).FPSM( JP) 

600  CONTINUE 

00  610  J=l . 10 
ETHTsETHT+ETH I  J) 

EPST=EPS  T+EP A  <  J ) 

ETMMT«ETHMT+ETHM( J) 

EPSMTaEPSMT+EPSM ( J ) 

610  CONTINUE 

IS  FORMAT  (5X . 33HTHE  FOLLOW  I NG  IS  FOR  TE  INCIO ENCE.Z) 

WRITE  (6.18) 

WRITE  (6.11)  ETHT.EPST 

11  FORMAT  ( lOX, lOHE-THETA  *  .2F12.6./. 10X.8HE-PSI  ■  .2F12.6) 

F)«4.*< J. 14159**5)* ( (0«W* A) 4<2I4<<CA8S<  I .-ER) )«2> 

S!GTH»F1#< (CABS(ETHT) )<<2) 

S I GPS*F  1  <  <  <  CABS  (EPS T  )  )«2) 

S IGTHD3 10. *AlOG1 0 ( S IGTH/REF ) 

S I GPSD* 10.*AL0Gl0(SI GPS/REF ) 

WRITE  (6.12) 

12  FORMAT  (/.8X.42HTE  SAME  POL.  BACKSCATTER 1 NG  CROSS  SECTION.) 

WRITE  (6.131  S I GPS . S I GPSD 

13  FORMAT  •( 15X.8MS1GMA  «  .F] 4.8 « 1 5X « 1 2MS1 GMA ( DB )  *  .F14.0) 

WRITE  (6.14) 

14  FORMAT  (/.SX.42HTE  OPP.  POL.  BACKSCATTER I NG  CR0S5  SECTION.) 

WRITE  (6«13)  SIGTH.SIGTHO 

19  FORMAT  (SX.33HTHE  FOLLOWING  IS  FOR  TM  INCIDENCE./) 

WRITE  (6.19) 

WRITE  (6.11)  ETHMT.EPSMT 
S IGTH*F 1 4 ( ( CABS (ETHMT ) )4#2 ) 

S 1 6PS*F 1 4 ( <  CABS (EPSMT ) )442 ) 

SIGTMD»l0.4ALOGlO(S IGTH/REF ) 

S I GPSD* 1 0 . 4 ALOG l 0 ( S I GPS/REF ) 

20  FORMAT (/.SX.42HTM  SAME  POL.  BACKSCATTERING  CROSS  SECTION.) 

WRITE  (6.20) 

WRITE  (6.13)  SIGTH.SIGTHO 

21  FORMAT  (/.SX.A2HTM  OPP.  POL.  BACKSCATTER I NG  CROSS  SECTION.) 

WRITE  (6.21) 

WRITE  (6.13)  SIGPS.SIGPSO 

WRITE  (6.18)  A.O.W.FREO.THO.ER.REF 

ETHT«(0. .0. ) 

EPSTXO..C.) 

ETHMT  a ( 0. .0 .  ) 

EPSMT® ( 0. .0. ) 

800  CONTINUE 

801  CONTINUE 
STOP 
END 

•IBFTC  CC  NOOECK 

SUBROUTINE  NRCELL(M.RC«PC.ZC .AL.NX1 .XM.YM.ZM.AA.BB.CC) 

NX!  IS  THE  MAXIMUM  PHI-MOOE  INDEX  PLUS  )• 

PC. PC  »Zr  DEFINE  THE  CENTER  OF  THE  NEAR  CELL 
XM.YM.ZM  ARE  The  MATH  POINT  COORDINATES. 

CELL  DIMENSIONS  ARE  R«AA .  PHI«BB.  L«CC. 

R|  THROUGH  P16  ARE  THE  VARIOUS  INTEGRALS. 

M  IS  THE  NUMBER  OF  SUBCFLLS  ALONG  EACH  EDGE  OF  THE  MAIN  CELL. 
M4M4M  IS  THE  TOTAL  NO.  OF  SUOCELLS.  M  IS  EVEN 
COMMON  /XXX/  XM ( 10) . ZM (10) 
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COMPLEX  XP.P.O 

COMPLEX  Pf  1  I  0 J  <P21 1 9  J  <P3<  J  OJ  <P4 1 1 0  J  <P51 10  J  <P61 1 0  I «P7 <  tO 1 <P8 ( lO 1 
COMPLEX  P9tlOi.P10<l0l.Pll<101.Pt2<t0l.P13<10l.P14<l<5j 
COMPLEX  P15< 101 .P16 ( 10) 

COMMON  PI ,P2.P0.P4. P5<P6<P7<P8 <P9<PI0<P1 i <P12.P13<PI4.P15.P16 

COMPLEX  XPP.XPO 

PEAL  L.LS 

00  00  N«1 .NX] 

t>l  <Nls(0<<0<  1 

02  INI" <0. <0<  I 

oj<Ni*n<  < o <  i 

P4<NJS<0<<0< 1 

P5<NJS<0<<0<J 

P6<NJS<0<<0<J 

P7<NJ*  <0<  <0<  J 

P6<N1S<0<<0<  I 

P91NJS<0<<0<J 

P10<NJS<0<<0<J  , 

Pit  <NJS<0<<0<  1 
P12<NJ«  <0 . <0. 1 
P13<NJS<0<<0<  J 
P141NJS<0<<0<  J 
Pl5(Nj*  < 0< <0<  J 
P161NJS<0<<0<  J 
50  CONTINUE 
TlsM 
T2sm/2 
AsAA/T J 

e«ee/Ti 

C*CC/T I 

RADS«RC-<T2-0.500 J*A 
PM I S»PC- <  T2-0 .500 1 *0 
l  5«-lT2-0.500l*C 
00  100  I s 1 <  M 

risi-i 

PMt*PHIS-*Tl*0 
C1«C0S1PM| 1 
4 1  *  S I N 1 PM 1  J 
OO  100  Jsl.M 
Tl*J-i 

RAOsPAOStTJ  *A 
00  100  Ks|  < M 

TtsK-t 
l *CS*T1*C 

X»XM- <PAO*COS 1 *L  J-L*SIN< AL J  J*COSIPMI  ) 

YsYN-<RAO»COS< AL J-l  »SIN< AL  J  J*SIN<PMl  J 

ZsZK-lZ^+LACOEtAL! 1 

P»SOPT <  x*X+Y* Y+Z*Z 1 

*3«R**3 

R5»p**5 

AQs«.2BJ19«P 

PRst-i .*AR«ARJ/R3 

O 1 s— AR/R3 

PsCMPLXtPR.Pl J 

0»«  1 3.<  -AB* AR  j  /R5 

0I«3<*AR/R3 

QsCMPLXlOR.Qt  J 

APsRAO»“0<jl  *L  J 

OV*AP*JJ*A  *0/6.283 19 

XPsCMPLX <  COS  I AR1<— S1N1ARJ  J 

XPPsXPSP 

xposxp*o 

00  I  00  N* 1 <NX1 
T2*N-1 


Fig.  24--(continued) . 
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CN«C.0S<T2*PH1 ) 

SN«SIN<T2*PH1 ) 

CNO»CN*OV 
SNO«SN*OV 
Oj INI *P 1 <N>+C1*CN0*XPP 
P2(NI«P2(N1+SI*SN0*XPP 
P3fN)»Pj (N J+X*X*C1 *CNO*XPO 
0*1N>«P4 (N)*X*X*S1*SN0*XPQ 
PS«N)*P3<N>+X*y*S1  *CND*XP<1 
06«n»«P6<N)+X*Y*C1*SND*XPQ 
O7(N»»P7<NM0r*Z*CND*XP0 
08  C N  )  »P8  «  N  I  ♦  Y*Y*S  1  * SND*XPQ 
09 1 N ) »P9  <  N  >  +Y  *Y*C 1 *CN0*XPQ 
OlO(Nl»PlO<N)+Y*Z*SNO*XPO 
Ol 1 (N)«P1 1 (NI+CN0*XPP 
O  12  <N  I  «P1  2  <N  > +Z*X*C 1 *CN0*XP0 
OJ3<N)«P13<N)+Z*X*S1*SNO*XPQ 
o  1 4  ( N  >  *P  1  4  < N )  «-Z*7*S  l  *CN0*XPQ 
P|S(N1»P15(N)  ••■Z*X*C  l  *SN0*XPQ 
0|6<N1»P16<N)  +Z*Z*CN’0*XPQ 
100  CONTINUE 
PETUPN 
END 

•DATA 

9  0.750  O.lOO  0.200  6.030 


Fig.  24— (continued). 


B.  Tube  Scattering  with  an  Axially  Incident  Plane  Wave 

The  flow  chant 

of  the  basic  operations 

is  shown  in  Fig.  24  and  the 

statement  listing  is 

are  defined  below. 

given  in  Fig.  25.  The 

input  and  calling  parameters 

The  data  cards 

needed  for  each  run  are 

as  follows: 

Card  Number  1 

Format  (2I5.7F10.5) 

Columns 

1 

Quantity 

Description 

i-5 

NNN 

NNN  =  the  maximum  ^-mode 

(right  adjusted) 

index  (NNN  =  2  for  on- 
axis  incidence) 

6-10 

MMM 

MMM  =  the  maximum  z-mode 

(right  adjusted) 

index. 

11-20 

A 

A  =  the  mean  ring  radius 
in  inches. 

21-30 

T 

T  =  the  tube  wall  thick¬ 
ness  in  inches. 

31-40 

TL 

TL  =  the  total  tube 

length  in  inches. 

41-50 

FR 

FR=  the  frequency  in 
gigahertz. 

51-60 

Re(ER) 

ER  «  the  complex  relative 
dielectric  constant  of 
the  tube  material . 

61-70 

Im(ER) 

Card  Number  2 

Format  (1115) 

Columns 

Quantity 

Description 

1-5 

JT 

JT  *  the  total  number  of 

(right  ad-*iisT.-d) 

rings  into  which  the 
shell  is  divided. 

76 

The  calling  parameters  for  the  RING  subroutines  are  defined  below: 
NX  s  the  maximum  $-mode  Index 
MX  *  the  maximum  z-mode  Index 
A  *  the  Inner  radius  of  the  tube  in  inches 
Z1  *  the  z-coordinate  of  the  center  of  the  particular  ring 
(In  waelengths) 

AL  *  0.0  (always) 

T  *  the  ring  thickness  In  wavelengths 

TL  *  the  total  tube  length  In  wavelengths 

DL  a  the  length  of  the  particular  ring  (in  wavelengths) 

The  calling  parameters  for  the  NRCELL  subroutine  are  defined 
below: 

M  *  the  number  of  subdlvldlons  of  the  large  cell  along  each  side. 
There  are  M3  total  subcells. 

RC»PC,ZC  a  the  p',$,  z  coordinates,  respectively,  of  the  center  of 
the  large  cell. 

AL  *  0.0  (always) 

NX1  *  the  maximum  #-mode  Index  plus  1. 

XM,YM,ZM  *  the  rectangular  coordinates  of  the  match  point  for 
a  given  calculation. 

AA,BB,CC  *  the  large  cell  dimensions  In  terms  of  the  polar  coordin¬ 
ates,  p,  4,  z  respectively. 

The  function  FL(Z)  gives  the  distance  along  the  shell  (tube)  as  a 
function  of  the  z-coordinate  of  the  point  on  the  shell.  This  generalized 
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function  approach  Is  used  In  anticipation  of  the  general  case  when  the 
arc  length  along  the  shell  will  be  a  more  complicated  function  rf  z  such 
as  a  polynomial  with  several  terms. 


Fig.  25— Flow  chart  of  the  tube  scattering 
computer  program. 
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RING  SUBROUTINE 


Fig.  25— (continued) . 
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Fig.  25--( continued) . 
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van  ooeren,  r.  e.  jobn  fgb230  a 

•EXECUTE  (8 JOS 

•IB JOB 

•18FTC  MAIN  NOOECK 

1  FORMAT  (215.7F10.5) 

22  FORMAT  (5X.AMJ  »  .19) 

2  FORMAT  (12F10.4./) 

3  FORMAT  (SX.6E 1 6.  7.  /  .  5X .  6E  1  6 .7  ./.5X  •  6E I  6.7. // ) 

4  FORMAT  (5X.6HMXX  ■  .19) 

5  FORMAT  I 5X.6E16.7./.10X.6EI6.7./) 

6  FORMAT (5X. J 7MER. T.L . A. WAV.FR  « 7F I  1 • 6 . / • 5X «8HS I GM A  •  .F16.9) 

7  FORMAT  1 2X. 10HATAUI0. I )«.2F10.S.1|M  BPH1 (0.1>«.2F10.5.1|H  ATAUI I • I  I 
2®.2F|0.5.1 1H  BPH 1(1.1)*. 2F 10.5) 

8  FORMAT ( 9X .2 1 9 ) 

9  FORMAT  I9X.9HJT  ■  .15./) 

10  FORMAT  (1115) 

11  FORMAT  (5X.I4HTHE  A.S  FOLLOW) 

12  FORMAT  (5X.14HTHE  B.S  FOLLOW) 

13  FORMAT  (5X.28MTME  AUGMENTED  MATRIX  FOLLOWS) 

14  FORMAT  (5X.28HTHE  AUXILIARY  MATRIX  FOLLOWS) 

COMPLEX  FTHS.FS 

COMPLEX  A | (6.3.2) • A2(6< 3.2 ) . A3( 6.3.2  > . A4 (6.3*2 )• 45(60.2) 

COMPLEX  A6(6. 3.2). A7(6. 3.3). AS (6.3.2 ).A9(6.3.2) 

COMPLEX  B1 (6.3.2) .82(6. 3.2) .83(6.3.2 ) .64(6.3.2) >B9 (6. 3.2) 

COMPLEX  86(6.3.2) .87(6. 3.2  > .88(6.3.2 ) .89(6. 3.2) 

COMPLEX  AA1 (6.3.2 X AA2( 6.3.2 )<AA3(6. 3.2). AAA (6.3.2 >»AA5 (6.3.2) 
COMPLEX  AA6(6.3.2) . AA7(6.3.2 ) . AA8 (6. 3. 2 ) . AA9 (6. 3.2 > 

COMPLEX  881 (6. 3 .2 )• 882 (6. 3.2). 863(6.3. 2 >.B84 (6. 3.2 >.885 (6. 3.2) 
COMPLEX  886(6.3.2) .887(6.3.2 ). 868 (6.3.2 ).BB9I6. 3. 2) 

COMMON  /AAA/  AA| , AA2.A43.AA4.AA5.AA6.AA7.AA8. AA9.BB1 .682.883.884.8 
285 . BB6 . 887 . BB8 . 889 
COMPLEX  C(72,72).E.FF.FTH.FPS.ER 
COMPLEX  ARG.ARG2.EX 
COMMON  /XXX/  XM(IO) .ZM( 10) 

COMMON  /JWW/  JWW 
DIMENSION  JZ(10) 

INTEGER  OOEV 

READ  (5.1)  NNN.MMM.A.T.TL'.FR.ER 
NNN  is  the  MAXIMUM  PH I -MODE  INDEX. 

MMM  IS  THE  MAXIMUM  L-MOOE  INDEX. 

A  IS"  THE  MEAN  RING  RADIUS  (  OR  TUBE  RADIUS). 

A.T.1L.  ARE  IN  INCHES.  FR  IS  IN  GIGAHERTZ. 

ER  IS  THE  RELATIVE  DIELECTRIC  CONSTANT  ICOMPLEX  IN  GENERAL). 

DO  IS  THE  APPROX.  LENGTH  IN  WAV.  OF  EACH  ELEMENTAL  RING. 

NX>*4NN+1 
MX®MMM+ 1 
MMX®2*MX-| 

c  MMX  IS  the  number  of  matching  RINGS  OVER  THE  Tube. 
mi«mmm 

READ  (5.10)  JT . ( JZ ( I ) . I « 1 . MMX ) 

C  JT  IS  AN  ODD  INTEGER  (THE  TOTAL  NUMBER  OF  RINGS  FOR  INTEGRATION). 

WAV®1 I .803/FR 
•■A/WAV 
T®t/WAV 
TL*TL/WAV 
T I  ■  JT 
0L«TL/T1 
00  41  J® I .MMX 
T1®JZ(J)-I 
ZM< J)®TI*DL40L/2. 

XM I J )  ■  A 
41  CONTINUE 

WRtTE  (6. I)  MMX. JT. (ZM< J) . J®l .MMX) 

WRITE  (6*1)  MMX. JT. (XM( J) «J®1 .MMX) 


Fig.  26--Tube  scattering  computer 
program  statement  listing. 


A L*0.0  82 

N«NX 

00  82  J*1 .MMX 

00  82  M*| .MX 

At  < J.M.N>«(0. .0. ) 

A2( J.M.N)*(0. ,0. t 
A3( J«M.N>*(0. .0. ) 

A4( J.M»N)*<0, .0. t 
AS( J.M«N>«(0. .0. ) 

A6< J.M.N) * ( t>. ,0. > 

AT IJ)MiN)XOi  .0.  ) 

A8 UiWiNIXOi  .0.  I 
A9< J.M.N)* (0. .0. > 

81 C J.M.N)*(0. ,C. ) 

S2(J«M«N>*(0. .0. ) 

B3« J.M.N) * (0. .0* I 
B«( J.M.N) *  <0. ,0. ) 

BSI J.M.N>*(0. .0. ) 

B6( J.M.N ) *  <0. .0. ) 

BT< J.M.N)*(0. ,C» 1 
BSI J.M.N) * (b. .0. ) 

B9( J.M.N) * (O. .0. ) 

82  CONTINUE 

C  THE  SYMMETRY  OF  THE  TUBE  IS  USED  SO  I NTEGRAT I  ON  IS  NEFOEO  OVER**** 

C  ****  ONLY  ONE  HALF  OF  THE  TUBE 

JT2*(  JT-fl  )/Z 
AI-A-T/2. 

OO  SO  J1«1.JT2 
JKWJI 
T 1  *JI 

2l«(Tl*0L)-0L/2. 

N*NX- 1 

CALL  RING  (N.MI.At.Zt.AL.T.TL.OL) 

N*NX 

OO  88  J«1 .MMX 
00  88  M«1 .MX 

At  I J.M.N >»A1 1 J.M.N)*AAl (J.M.N) 

A2( J.M.N) *A2( J.M.N) +  AA2I J.M.N) 

A3< J.M.N >«A3( J.M.N )+AA3( J.M.N) 

A4( J.M.N) *A4( J.M.N )+AA4( J.M.N) 

AS( J.M.N) *AS( J.M.N )+AAS( J.M. N) 

A6(J«M.N>*A6<J.M.N)+AA6 (J.M.N) 

ATI J.M.N) * AT < J.M.N )*AA7 (J.M.N) 

AB( J.M.N )«A6( J.M.N )+AAB( J.M.N) 

A9( J.M.N) *A9( J.M.N I +AA9( J.M.N) 

B1  I  J.M.N)  *81  <  J..M.N  )+B81  (J.M.N) 

B2IJ.M.N) *82 (J.M.N) +BB2 (J.M.N) 

83 ( J.M.N) *83 (J.M.N ) *3B3l J.M.N) 

B4 ( J.M.N) *34 (J.M.N ) TBB4 (J.M.N) 

8S( J.M.N) *B5( J.M.N) .BBS (J.M.N) 

BA  I J.M.N) *B6( J.M.N) .8961 J.M.N) 

BT( J.M.N) «B7( J.M.N )+BB7( J.M.N) 

BSI  J.M.N)  *BB(  J.M.N)  -.000  (J.M.N) 

89 ( J.M. N) *89 (J.M.N ).BS9( J.M.N) 

88  CONTINUE 
SO  CONTINUE 

00  60  J*t .MMX 
DO  60  M*1 .MX 
WRITE  (6.8)  J.M 
WRITE  (6.11) 

WRITE  (6.8)  41 ( J.M.N). A2( J.M.N). A3( J.M.N) >A4 (J.M.N) «A8( J.M. Nl «A 

26 (J.M.N)«A7IJ.M.N).A8(J.M*N).A9(J.M.N> 

WRITE  (6.12) 

WRITE  (6.8)  81 (J.M.N).B2( J.M.N) .031 J.M. N) .04 (J.M.N ). BSI J.M.N). 8A(J 
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2.M.N).B7(J.M'N).B8<J.M.N).B9(J.M,N) 

60  CONTINUE 

e«er-i. 

C  NEXT  STEP  IS  CALCUL AT I ON  OP  THE  COEFFICIENT  MATRIX. 

00  70  Jt>l.MMX 
00  70  Ml  at .MX 

THE  FOLLOWING  SETS  UP  THE  COEFFICIENT  MATRIX. 

THE  MATRIX  IS  FILLEO  IN  FROM  LEFT  TO  RIGHT. 

THE  SMALLEST  VALUE  OF  M  IS  TO  THE  LEFT « M  INCREASES  TO  THE 
M«M1-1 

T  I  *M 
J«J1-1 

CM*C0S(6.2831 9»T1*ZM( Jl )) 

SM«SIN(6.283|9*TI*ZM< J1 ) 1 
I S«6#J 

IF  (MI .60.1 »  KS«0 

IF  (Ml.GE.2l  XS«12*M-6 

I t«IS*l 

I2«IS*2 

1 3* I S43 

I4«IS+4 

I S«t STS 

16*1 5+6 

X1«XS+1  . 

X2*XS+2 

X3»XS+3 

X4*KS+4 

XS*XS+S 

K6*KS+6 

X7*XS+7 

K8*XS+8 

X9*XS+9 

XlOaXS+lO 

Xll*KS+ll 

X 1 2*XS+ 1 2 

C(ll »Xl I«CM-E*AI ( Jl .Ml .N1 
C(I|.<2)«(0.,0.» 

C((|.K3)a-E«A3<Ul.M| ,NI 
C(l!«K4)a<0.,0. t 
C < 1 1 «XS1*-E*A2( JI .MI ,N) 

cm  .<6 >*  to.  ,0.  > 

IF  (Ml .EO.l »  GO  TO  SOO 
C(Il »K7)«SM-E*B1 ( Jl .Ml .N1 
C(lt »X8l* (0. ,0. * 
cm  »X9 ) «-E*B3 ( Jl  .Ml  ,N1 
C (11 .XI 0). (0..0. ) 

Cdt.Kl  |  )«-6*B2(  Jl  .Ml  ,N) 

Cm.K12)«(0..0.) 

900  CONTINUE 

C  <  I2.KI  IX0..0.1 
C ( I2.X2 1»C(II«<I I 
C(I2<K3|X0..9. ) 

C ( I2.K4 )»-C ( I  I .X3) 

C ( 1 2.XS 1X0..0.) 

Cl  I2.K6)*C( 1 1 .X5> 

IF  (Ml  .EO.l  )  GO  TO  501 
C ( 1 2«X7 ) X 0 . • 0.  ) 

C  < I2.K8l.CC  I  1 «X7 1 
C(I2.K9).I0..0.1 
C  (  I  2.X  1  0  )  •— C  (  1 1  . X9  1 
C ( I 2» XI 1 )a(0.«0. I 
C(l2.KI2>>Cm,XI  1  ) 

SOI  CONTINUE 

C ( I 3«XI  >«E*A4( Jl .Ml «  N 1 
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CIl3.l2)«(0..O.) 

C(l3«<3)«CM-e**6( J] .Ml .  n  ) 

C  I  I  3. <4 >■ (0* . O. ) 

C  < 1 3.KS1 «E*A3( Jl , Ml  , N ) 
e< l3*<6 )« (0. ,0. ) 

IE  (MI ,£0. 1 1  GO  TO  502 
Clt3«K7)«E*B4( Jl.Ml  ,N) 

C(l3.K8)»(D.,0.) 

C  (  I3«K9)»SM-E  #86  (  Jl  ,M1  ,N1 
C  < ! 3.K1 0>«(0, ,0. ) 

CU3.K1  1  )aE*ei ( Jl  .Ml  «  N  1 
Cll3«Kl2)»(0,,0, i 

502  CONTINUE 

C  I  1  4 «K  1  )a(0.,0.1 
C  tl4.K2l-.-Ct 13.K J 1 
Cll4.K31«(0.,0. ) 

Ctl4<K4)-C(t3.K3) 

C  IUiKS)«(0.,0,  ) 

C<  I  4  «K6  1  »-C  (  13. *51 
IE  (Ml .E0.1 )  GO  TO  503 
Ctl4.K7)»(0.,0.) 

C(l4«K8)»-C(I3«K7) 

Cl I4.K91-I0..0.) 

Cl  U.K|OI«C<  tO.KO  1 

c  <  I4.ki  n«(o»  «o»  i 

C  1 1 4  «K  1 2  1 »— c ( 1 3.KI 1  1 

503  CONTINUE 

C  < 1S.K1  )«-E*A7( Jl . M] .N1 
C 1 I  5«K2 1 « ( 0 • , 0. 1 
C ( I5.K3 )»-E*A9( Jl .Ml .Nl 
C I  I5.K4 1« (0. . 0. ) 

C I l5«<51«CM-E*Aft( Jl .Ml  ,n> 

CIIS.K6I«(0.,0.) 

IE  (M1.EQ.1I  GO  TO  504 
CII3.K7)«-E*87f Jl .Ml .N) 

C I  I  3  .K8 ) ■ ( 0 . • 0 . ) 

C I I5.K9)«-E*89( Jl .M|  .N) 

C I IS.Kt Ol-IO. ,0,  ) 

C I  1 5 .< I  I  |«SM-E*88( Jl .Ml ,N1 
C(l5.<12)«10,,0. ) 

504  CONTINUE 

C( Ift.Kl  1-I0..0. > 

C ( I6.K2 l-CI I5.KI 1 
C<I6.K3)«(0.,0, ) 

Cl 1  6  *K4 1  --C ( I3.K31 
Clt6.K5>-(0.,0.) 

C(I6.K6)«C( I5.K3I 
IEt  M1.E0.11  GO  TO  505 
Cl 16«K7)-(0.,0. ) 

CII6.K8l-CM3.K7> 

Clt6.K9>-(0..0,> 

CII6.K101--CI 15. <9 1 
C I I6.KI 1  I ■ (0.  ,0.  1 
CII6.KI2laCII3.IClI  I 

505  CONTINUE 
70  CONTINUE 

MX  I  «  ( 6*MMX 1 ♦ 1 

DOWN  TO  80  is  CALC •  THE  LAST  COLUMN  OK  THE  AUGMENTED 
00  80  J| -1 .MMX 
J4JI-I 
I S-6*J 

AB03«6.28319*2M( Jl  ) 

ABGaCMPuX (0. . A9G3 I 
I1«IS+1 
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12*IS+2 
I3-IS+3 
!4*IS+4 
1 5s 15  +  5 
1 6* ] S+6 

Cell «MX1 )=CEXP(-ARG I 
C ( I  2 • MX 1  >  =  10. .0. > 

Cl  I3«MX1  )=>-CEXP(-ARG) 
cii4.Mxn*io..o.) 

C I  I 5.MX1  ) * (0.  .0.  ) 

C I I6.MX1 ) aio. .0. j 

1 F  (J.EQ.l)  WRITE  (6.1)  J . J • ARG.C (1 .MX  1  ) ,C 1 3. MX1) 

80  CONTINUE 
JJJ«6*MMX 
JJJ1»JJJ+1 
WRITE  (6.13) 

00  82  J*l .  JJJ 
WRITE  (6.8)  J 

WRITE  (6.5)  (C(J.M).M»1.JJJ1 ) 

82  CONTINUE  ! 

THE  COEFFIECENT  MATRIC  IS  6*MMX  BV  6*MMX. 

THE  AUGMENTED  MATRIX  IS  6*MMX  BV  (6*MMX+1 1 

THE  FOLLOWING  IS  THE  CROUT  REDUCTION  OF  THE  COEFFICIENT  MATRIX. 

MXX*6*MMX 

NN«MXX+1 

DO  1  IB  L  =  1 . MXX 

l  LL*L-1 

DO  118  1C-L.MXX 
1 I=1C+1 

IF  (LLL.EQ.O)  GO  TO  1171 
DO  117  X* J .LLL 

C(IC.L)«C(IC.L)-C(  IC.O*C(K,L) 

I  17  CIL.I1  >=C(L.I1)-C(L.K)«C(K.1I  I 
1171  CONTINUE 

118  C(L. I  I  >»C(L. 1  I I/CIL.L) 

WRITE  (6.14) 

WRITE  (6.4)  MXX 

00  120  J*l .MXX 

WRITE  (6.22)  J 

WRITE  (6.2)  (CIJ.I I. 1*1. MXX) 

120  CONTINUE 

IF  IMXX.Eff.l)  GO  TO  124 
DO  123  (.*2* MXX 
IC-NN-L 
I I=IC+1 

DO  122  K«I I .MXX 

122  C ( IC*NN)«C( IC.NN)-C ( IC.X)*C (K.NN) 

123  CONTINUE 

124  DO  126  J»1 .MXX 
WRITE  (6.3)  C(J.NN) 

126  CONTINUE 

THE  FOLLOWING  COMMENTS  APPLT  FOLLOWING  THE  CROUT  REDUCTION. 

C(l «NN)»A(0.1 )  SUPERITAUI 
C (2.NN ) *B (0  *  1  )  SUPER  (TAU1 
C ( 3.NN ) *B(0. | ) SUPER (PHI ) 

C(4.NN)«A(0.1  ) SUPER (PH  I  ) 

C(5.NN)*A(0.1 )  SUPER  <L) 

C (6.NN)»B(0. 1 )  SUPER  (L) 

C (M.NN)«AS(0. 1 )  SUPER  (TAU) 

C ie.NN)»BSIO, 1 )  SUPER  (TAUI 
AND  SO  ON 

NN  IS  THE  INDEX  OF  THE  LAST  COLUMN  OF  THE  REDUCED  MATRIX. 

WRITE  (6.7)  C ( 1 .NN) .C ( 3.NN1.C (7.NN) .CI9.NN) 

C  the  REST  IS  calc.  OF  THE  SCATTERING  CROSS-SECTION. 
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FTH*<0.,O.I 
FPS« (0. .0. ) 

FTHS  « (0. .0. ) 

JTTTT«2*JT 

0L«0L/2. 

00  201  M«1 <MX 
FF*(0..0. ) 

FS«<0..0. ) 

MM.M-1 

T2«M-1 

IP  (M.EO. I I  MS*0 

IP  CM.GE.2)  »S«12*MM-6 

MS1«MS+1 

MS2»MS*2 

MS3»MS+J 

MS4»MS*4 

MST«MS+T 

MS9«MS*9 

00  200  J« 1 «  JTTTT 

T  I  ■  J 

Z«(T1*0L I-0L/2. 

ARG1-6. 28319*Z 
ARG2«CMPLX<0. .-ARGI I 
EX«CEXP<ARG2) 

CM«C0S(2.»T2*3. 14159*21 
SM«SIN12.*T2*3.I4199»Z) 

FF«FF+CM*EX 

PS«FS*Sv|*EX 

200  CONTINUE 

rTH.FTH-MCIMSI .NNI-C1MS3.NN) I *PF 
IP  (M.EO.l )  GO  TO  201 1 
pTMS«FThS  ♦1C<MS7.NN)-C(MS9.NN>  >*PS 
*011  CONTINUE 

•BITE  16.21  PP .FS.FTM.PTH5 

201  CONTINUE 

C  NN  IS  THE  INOEX  OF  THE  LAST  COLUMN  OF  THE  REDUCED  MATRIX. 

CO  I «4 . *  1 3. 1 4  I  99**8 1  * ( I T*OL*A  >  **2  > 

C02* <  CABS  1 1 .-ER 1 )**2 

SIG.C01 *C02*< (CABS1FTH+FTHS I  1**2) 

WRITE  I6.6>  ER.T.U.A.WAV.FR.SIG 

202  CONTINUE 
STOP 
ENO 

•*iBPrc  oecki  '  nooeck 

SUBROUTINE  RING  (NX . Ml , A , Z 1 . AL .T . TL . DC ) 

C  ALL  DIMENSION  ARE  IN  WAVELENGTHS. 

C  NX  IS  THE  MAXIMUM  PHI  MODS  INDEX. 

C  Ml  IS  THE  MAXIMUM  l.-MOOE  INDEX. 

C  21  IS  THE  Z— COORDINATE  OF  THE  CENTER  OP  THE  INNER  SURFACE  . 

C  AL«0.0  IN  ALL  CASES. 

C  AL  IS  THE  VAH«  OF  ALPHA  AT  THE  RING  BEING  INTEGRATED  OVER  HERE. 

C  A  IS  THE  INNER  SURFACF  RADIUS  OF  THE  RING  BEING  INTEGRATED  OVER. 

C  T  IS  THE  PING  THICKNESS  NORMAL  TO  THE  SHELL  SURFACE. 

C  OL  IS  THE  INCREMENTAL  ARC  LENGTH  ALONS  THE  SHELL. 

BEAL  L.LJ 

COMMON  /XXX/  XOl |0| .201 10) 

C  X0.20  ARE  THE  MATCHING  POINT  COORDINATES. 

COMMON  /JWW/  JWW 
MX*M I ♦ 1 
MI»**MX-I 

C  Ml  IS  THE  NUMBER  OF  MATCHING  BINGS  OVER  THE  TUBE. 

NX 1 bNX.1 

OIMENSION  ST(9 | «  SPI 9 1 . SL (9  I 

COMPLEX  XP.P.O.PI (3) .P213) .P313).P4(3> .P913) .P613) .PTt3 1 .oo (31 
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COMPLEX  P9t31.PlOI3l.Pllt31.Pl2(3!.P|3(31.P!4<31.Pl5(31.P16(31 
COMMON  Pl.P2.P3.P4.P5.P6.=7.P3.P9«Pl0.PII.PI2.P|3«P14.Pt5.P16 
COMPLEX  V](5.2l.V2t5.21.V3tS.21.V4«5.21.VS(S.2).V6<5.2I.V7(S.2l 
COMPLEX  va(5.2).V9(5.21«V10(5i2I.Vll<5.2t.V12(S.2>.V13(5.2) 

COMPLEX  VI4(5,2),VI5(‘5.2).V!6(S.2) 

COMPLEX  XPP.XPO 

COMPLEX  41 (6. 3.2) . 42  <6. 3.2 ) . 43(6. 3.2 1 • 44 (6. 3 .2 I . 45(6.3*2) 

COMPLEX  46 <6. 3, 2). 47(6. 3*2). 48 (A. 3. 2). 49 (6. 3. 2) 

COMPLEX  8)  (6. 3. 2>. 82 <6. 3. 2  J .83(6. 3. 2  >.04 (6. 3. 2). 85 (6. 3, 2) 

COMPLEX  86 <6. 3. 2). 87(6. 3. 2). 88(6.3. 2). 89 (6. 3. 2) 

COMMON  /44A/  41 .42.43.44.48. 46.47. 48. 49.B1 .B2« 83. B4. R6.B6. 87.BQ.B9 
FIRST  SUBSCRIPT  IN  A1  ETC.  IDENTIFIES  THE  MATCHING  PING. 

SECOND  SUBSCRPT  IS  THE  L-MOOE  INDEX  PLUS  1 

THE  THIRD  SUBSCRIPT  IS  THE  PH  I -MODE  INDEX  PLUS  I. 

00  SO  N*  1 .NX! 

OO  SO  1*1 .MI 
VI ( I «N)  *  <0. .0. ) 

V2(I.N)  * ( 3. « 0, > 

V3tI.Nl  *tO.,0.1 
V4<(.NI  *(0. .0.) 

vsci.N)  *<o..o.t 

V6II.N)  *(!).. 0.1 
V7ll.Nl  *(0..0.l 
VS(I.N)  *10. .0.) 

V9(I.NI  *(0..0.1 
VlO< 1 .N»*(0.,0, 1 
VI  1  (  t «N1* (O. ,0. 1 
V12< I .Nl* (0. .3.  1 
V13»I.N)*(0..0.) 

V  I  4 ( I .Nl* (0. .0. 1 
VlSf l.N|*<0..0. 1 
V16( I .Nl* (0. .0. 1 


so 

CONTINUE 
OO  SS  J* 

1  «M] 

DO  S3  M* 

1  *MX 

00  SS  N* 

I  9  NX  1 

At ( J.M.N 

)*  (0. 

.0. 

42( J.M.N 

)*  (0. 

.0. 

431 J.M.N 

)«(0# 

.0. 

A4( J.M.N 

>«  (0* 

.0. 

4S (J.M.N 

)*<0* 

.0. 

46 ( J.M.N 

)  *  (0* 

.0. 

A7< J.M.N 

)  *  ( 0  • 

.0. 

A8t J.M.N 

)*<0, 

.0. 

49* J.M.N 

>*(0. 

.0. 

SI (J.M.N 

1  *  ( 0* 

.0. 

B2< J.M.N 

)*  (0, 

•  0. 

B3( J.M.N 

)  •  <  0. 

.0. 

SAC J.M.N 

)  *  (  0  # 

.0. 

B3< J.M.N 

)*  (0. 

.0. 

B6( J.M.N 

I»  <0# 

.0. 

B7t J.M.N 

>■<0. 

.0. 

BBC J.M.N 

)i(0« 

.0. 

B9( J.M.N 

)  *  (0. 

.0. 

ss 

CONTINUE 

00  tOO  J«t .Ml 

J  HERE  RUNS  OVER  THE  NUMBER  OF  MATCHING  RINGS  ON  THE  BOOT  (SHELL). 
ALPHA  MEASURED  FROM  TANGENT  TO  Z-AXIS.  CCW  IS  POSITIVE.  CW  IS  NEG. 
RC« 1 4/C05  <  Al 1 14T/2. 

AC»RC*C0SIAL1 
ZC«ZtT «T/2. 1»SIN(4L ) 

XC*AC 

N0«3.I AI59*4C/0.0I 
IE  (NO.I.T.30)  NO*  30 
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TT|»NO 

0PHI«3.|4159/TT1 
OV*T*OL • AC*DOH | /6. 283 1 9 

RTEST«SORT< (XC-XOI J) ) »*2+ < ZC-ZO 1 J) >**2> 

00l«l.l*r>L 
002*2 • 1 *0L 

IF  (RTEST.LT. 002)  K0*11 
IF  <9TEST.LT. 001)  K0*16 
tF  <9TEST.LT. 0.001)  K0*22 
IF  (RTEST.GT.D02)  <0*1 

C  00WN  TO  200  IS  CALCULATED  THE  CONTRIB.  OF  THE  FAR  PARTS  OF  RING. 

N01 *NOT | 

00  200  K*KO.NO| 

TT|*K-1 

PM«(0PHI/2, )tTTI*OPH! 

C*COS<PH) 

S*SIN(PHI 
X*XO< J)-AC*C 
V*-AC*S 
2*201 JI-ZC 

IF  (Z.LT.O.OI)  w1)T«0.666667#<  (OL/2.  )**3> 

IF  (Z.GT .0.0) )  VTT«Z*Z 
R»S0RTIX»*2*V**2+Z**2) 

AP*6.28J1 9*R 

ARS*AR*#2 

R3«R**3 

pb*r»«s 

OR- < - 1 . TARS  1 /R3 
O  J *-AR/R3 
P*CMPLX(PR.Pl 1 
OR* (3.— ARS1/R5 
0I«3.*AR/RS 
0«CMPLX<QR.Q1 ) 

XP*CMPLX < COS < AR ) , -s 1 N < AR ) 1 

XPP«XP*P 

XPQ*XP*0 

00  200  N*l • NX  1 

TT2*N— 1 

CN*C0S<TT2*PH) 

SN*SIN(TT2*PH) 

VI  I J.Nl-VI ( J.NM-CN*C*OV*XPP 
V2  ( J.N )  *V2 <  J « N  1  -fSN*S*OV*XPP 
V3<  J;N)*V3<  J.NI+X*-X*CN*C»DV*XPa 
V4< J.N)«V4 ( J.NI+X*X*SN*S*DV*XPQ 
V3< J.NI*VS< J.N>*X*Y*CN*S*DV*XP0 
V6< J.N>*V«( J.N)+X*Y*SN*C»DV»XPO 
V7( J.N>*V7( J.N)+X*Z*CN*OV«XPQ 
V8 I J • N  >  * VB  <  J • N  >  Y*Y*SN*S*OV*XPO 
V9I J.N)*V9< J.N)+Y*Y*CN*C*OV*XPQ 
VI 01 J.N)*V)0( J,N)TY*Z»SN*OV*XPO 
VI  I  I J.N ) *V 1  I  ( J.N)TCN*OV*XPP 
VI2I J.N)*Vl2( J.N)TZ*X*CN*C*OV*XPO 
VI 31 J.N)*V13< J.N)+Z*X*SN*S*OV*XPO 
VI 4  I J.N1*V14< J.N) +Z*Y*CN*S*OV*XPQ 
VIS  I J.N)*V1S( J.N)*Z*Y*SN*C*OV*XPQ 
V16(J*N)»V16(  J.N) ♦VTT*CN*OV*XPQ 
200  CONTINUE 

00  201  N*l*NXl 

1  FORMAT  (//.5X.4IS) 

2  FORMAT  (SX.8F1I  .6./.SX.8F1  I  .6./.5X.SFI  1.6./.3X.8FI  | 

WRITE  <6*1 >  J.N 

WRITE  <6«2)  VI (J.N) ,V2< J. N ) . V3< J.N) , V4 < J.N ) . VS ( J.N > , VA ( J.N).V71J.N 
F).V8<J.N).V9<J.N).V|0<J.NI.VII<J.N).VI2<J.N).VI3IJ.N),VI4<J.N),v1S 

2< J.NI.V16I J.N) 


Fig.  26- -(continued) 
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c 

c 


20  I  CONTINUE 

IF  IKO.EO. It  GO  TO  71 
IF  IKO.EO. Ill  N1 *2 
IF  IKO.EO. 16)  Nl*3 
IF  IKO.EO. 22)  Nl«4 
NOTE  THAT  KO* 1  OR  Nl«2,3«4. 

IF  KO*|.  THE  MATCH  POINT  IS  FAR.  NRCELL  SUBROUTINE  IS  NOT  NEEDED. 
DIMENSION  MZI4) 

IF  INI .ME. 2)  GO  TO  61 

MZtl>*4 

MZI2>*4 

GO  TO  6J 

61  IF  1N1.NE.3)  GO  TO  62 
MZI 1 )«l O 

MZ 1 2 1*4 
MZI3)*4 
GO  TO  6J 

62  MZI) ) *20 
MZI2  )* I  0 
MZI3)*4 

MZ I  4 ) *4  I 

63  CONTINUE 

DO  70  1*1 .NI 

TT3»t-l 

OPH«3.*DPHI 

IF  INI.NE.4)  GO  TO  63 

OH I *OPH I +DPH/2 . +TT3*DPH 


GO  TO  66 

63  OHI*0PH/2.+TT3*DPH 

66  CONTINUE 

XC*AC*COSIPH| ) 
VC»AC*SINIPHI > 
ZC*ZC 
XM-XOI J) 

YM*0.0 

ZM-ZOIJ) 

•  A*T 
BB*OPH 


CC*OL 
MM»MZ I  I  ) 

CALL  NOVELL! MM. AC. PH I .ZC.AL.NX1 .YM.YM.ZM, AA.BB.CC ) 
DO  70  N«l .NX t 

N)+PI IN) 

N )  +P2  I N ) 

N  )+P3 1 N ) 

N ) +P4 IN) 

N)TP3 IN) 

NI+P6IN) 

N)*P7IN) 


VI  I J . N ) *V 1 <J 
V2I J.N)*V2! J 
V3< J.N)«V3I J 
V4I J.N)*V4( J 
VS< J.N)*V3I J 
V6I J.N)*V6I J 
V7IU.N)*V7I J 


V8I J.N)»V8I J 
V9< J«N)*V9< J 
VlOlU.N)*VIO 
VII  I J.N ) *VI 1 
V12I J.N)*VI2 
V 1 3  I J.N ) »v 1 3 
VI4I J.N)*VI4 
V  I  3  I  J.  N  )  *v  1 3 
V 1 6  I J.N ) *VI 6 
WRITE  16.11 
WRITE  16.2) 
20 1 0 1 N )  . P I  |  IN 
WRITE  <6.2 > 
2) .V8< J.N) . V9 


Fig,  25  —  ( continued ) 


NI4P8IN) 

N )  +P9 1 N ) 

J.N) *P | 0  I N ) 

J.NI4P1 I  IN) 

J.N) 4PI 2  IN ) 

J.N )  4P|  3  IN ) 

J.NI4PI4IN) 

J.N)4P|3IN) 

J.NI4P16IN) 

.J.N 

Pi <N).P2<N),P3<NI»P4INI.P3IN).P6IN).P7IN),P8IN>.P9IN>, 
.PI2IN).P)3IN).P|4IN).P1SIN).PI6IN) 

V) IJ.N) .V2<J,N).V3IJ.N), V4 I J.N ) *  V3 I J. N ) • V6<  J.N). V7IJ.N 
J,N)«V10(J,N).V11<J,N),V12(J«N).VI3IJ.N),v14(J.N).V13 


Fig.  26— (continued) 
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2CJ. N).V16(J.N) 

70  CONTINUE 
7t  CONTINUE 

IF  (KO.LT.22)  GO  TO  60 

IE  K0«22.  The  matching  POINT  IS  IN  THIS  RING  ITSELF. 

DOWN  TO  80  IS  THE  SINGULAR  CELL  CONTR I SUT ION. 

A  1  •T/'2« 

81»AC»0PHI 
CI*0L/2. 

SO*SORT  ( A I **2+B 1 **2+C 1 **2 ) 

TT2a2./J.I4lS9 

*T  < JJ«-I .tTT2*(ATAN< (AI*C1  )/(Bi*50)  >tATAN( ( A 1*81 >/(Cl*SQ) >  > 

SP( J)a-| .TTT2*( ATANI (BI*AI )/<CI*SO> >+ATAN( (Cl *81 >/(Al *SO> > ) 

SL« J)*-l .+TT2«( ATANI <CI«B1 )/<A|*SOt >-.4TAN( (CI*AI )/'(Bl#SO) > > 

3  FORMAT  (5X.8F1 I .6) 

WRITE  16.1  I  J 

WRITE  (6«9  >  ST(J).SP(J).SL(J> 

GO  TO  83 
60  CONTINUE 
ST  <  J  J  *0. 0 
SPC J>«0.0 
SL  ( J )  ■  0  •  0 
66  CONTINUE 
N«2 

6  FORMAT  (3110) 

6  FORMAT  (8F10.6) 

PUNCH  3.JWW.J.N 

PUNCH  6. ST ( J ) * SP( J)  . SL (J).VS ( J. N ) • V2 IJ*N)*V3(J.N)*V4(J*N) . VS ( J.N ) » 
?V6< J.N).V7( J.N).V8( J.N) .V9( J.N) «V|0(J.N).V1 ](J.N).V12(J.N).VI3(J.N 
2).V14(J»N).V13(J.N).V16(J.N) 

C  JWW  IS  THE  RING  INDEX. J  THE  M.P.  INDEX. N  THE  PHI-MODE  INDEX. ** 

C«COS(ALP( J) ) 

S«SIN(AlPU)  1 
CA«COS(AL> 

sa«sin(al i 

300  CONTINUE 
l fFLIZI I 

c  fl (2 1  is  a  function  giving  l  As  a  function  of  z  . 

DO  90  Na 1 . NXt 
DO  90  Ma 1 .MX 
TT9»M-I 
6«TT9*3. 14139 
CM«C0S(2.*G*L) 

SM»SIN(2.*G*L) 

COMPLEX  TI .T2.T3.T4.TS.T6.T7.+8.T9.T10.TI 1.TI2.T13.T1A.T13.TI6 
Tl»VKJ.N) 

T2*V2 (J.N) 

T3«V3 ( J.N ) 

T4«V4(J.N) 

T9«V3(J.N> 

T6*V6( J.N) 

T7*vT ( J.N I 
T8«V8 ( J.N I 
T9*V9 ( J.N ) 

T10>VIO( J.N) 

TllaVIt (J.N) 

T12"V12( J.N) 

T 1 3»VI  3 ( J.N ) 

T I 4aVl 4 ( J.N ) 

TI9«VI3(J.N) 

T(6«VI6(J.N) 

A) ( J.M.N) »CM#  <  T I .T3-.T3TST { J  > ).A )  ( J.M.N) 

A2 ( J.M.N) »CM*T7*A2 ( J.M.N I 

A3( J.M.N )»CM*(-T2-TS+T6 I +A3( J.M.N) 


Fig.  26— (continued) 
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I 

i 


A 


uuuuuu 


*4 ( J.M.N! *CM» ( T24T6+T8 >.A« ( J.M.N) 

AS ( JtMtNI «CM*T | O+AS ( J.M.N) 

46 ( J.M.N)*CM*  (TI-T54T94SP  ( J)  !4A6( J.M.N! 

47 « J.M.N) *CM* (T1 2+T14 |4A7< J.M.N) 

AS( J.M.N) »CM* (T1 1 +T 1 6+SL ( J 1 )+A8( J.M.N! 

A9( J.M.N) «CM*  <-TI 3+  T| 5 I4A9( J* M.N1 
8  I  (  J.M.N)  «SM*  (T  I  4T34T54ST  ( J  )  >+Bl (J.M.N) 

Bit J.M.N! *SM»T74B2(  J.M.N) 

B3<  J.M.N)  «SM4  (-T2-T4*T6)4B3(  J.M.N) 

B4 I J.M.N) »SM« ( T24T6+T8 I +B4 ( J.M. N ) 

B5( J. M.N) «SM»TlO+B5< J.M.N) 

B6(J«M,N)*SM*(T1-TS4T94SP(J)  )+B6(J.M.NI 
87  C  J.M.N  >*SM*(T|2+T1«)  *<?7<  J.M.N) 

08 ( J.M.N) »SM* (T I 1+TI6+5LI j) >+B8( J.M.N) 

B9(J.M.NI«SM*<-T134T15)489( J.M.N) 

C  THE  FOLLOWING  USES  The  MIRROR  IMAGE  SYMMETRY  TO  CALC.  ALl.  THE  4.S 

IF  ( (4BS(2l-(TL/2. ) )).LT. O.OOI )  GO  TO  90 
JJ«M1-J4| 

t  J»TL-L 

CMJ«COS(2.*TT5*3. 14159*LJ) 

SMJ«SIN(2.»TT5»3. 14 I59*LJ> 

41  < JJ.M.N)«CMJ*(T|4T34T54ST( J)  1441  <  JJ.M.N) 

42 < JJ.M.N)»CMJ*(-T7>4A2{ J J.M.N I 

A3( JJ.M.N >»CMJ*( -T2-T44T6  >4A3( JJ.M.N) 

A4I JJ.M.N 1'CMJ* ( T24T64TS ) 4 A4 (JJ.M.N) 

45 < JJ.M.N )«CMJ»I-TI 0  >4A5( JJ.M.N) 

46 (JJ.M.N | *CM J#  < T l — T5+T94SP I J ) > 4  A6 (JJ.M.N! 

47( JJ.M.N )*CMJ*(-T I 2-T 14 )4A7< JJ.M.N! 

A8(JJ.M.N)'CMJ*(T|  1  4T  1 6+SL  (  J  1  )4A8(JJ.M.N) 

A9( JJ. M.N! «CMJ*(Tl3-T 151+A9I JJ.M.N! 

B I ( JJ.M.N)»SMJ*(T14T34T54ST( J) 1481 (JJ.M.N! 

B2( JJ.M.N )»SMJ*(-T 7! +02 (JJ.M.N! 

83 ( JJ.M.N | »SMJ*(-T2-T44T6>4B3 (JJ.M.N) 

B4 ( JJ . M .N 1 *SM J* ( T24  T64T8 1 484  <  JJ.M.N ) 

BS( JJ.M.N |*SMJ»(-Tl 0 >4B3( JJ.M.N I 

B6( JJ.M.N )*SMJ*(T1-T34T94SP( J 1 >486  ( JJ.M.N! 

87 ( JJ.M.N )»5MJ*(-T1 2-T 14 )4B7( JJ.M.N) 

B8( JJ.M.N )«SMJ*( T1 1 4T 1 6+SL ( J ! 1488 ( JJ.M.N! 

89 ( J J.M.N! »SMJ*(T13-T1 5) 4B9 ( JJ.M.N) 

90  CONTINUE 

100  CONTINUE 

8  FORMAT  I5X.37HTHE  *«S  (FROM  RING  SUBROUTINE)  FOLLOW ) 

9  FORMAT  I3X.3~HTHE  B.S  (FROM  RING  SUBROUTINE!  FOLLOW) 

10  FORMAT! 3X. 31 3.6E 16. 7./.5X . 6E 1 6.7./.5X . 6E 16.7./ ) 

DO  101  I«t .Ml 
DO  101  J* 1 . MX 
DO  101  K-1.NX1 
WRITE  (6.8! 

WRITE  I6.|0)  |,J.K.Al<I.J.K).A2(I.J.K!.A3(l«J«K>.A4(f.J.K).AS(|.J 
2K ) . 46 ( I «J.K),A7( I .J*K).A8( I .J.K ).A9(! ,  J.K  ) 

WRITE  (6.9) 

WRITE  (6.1°)  I .J.K. 81 ( I .J.K 1.82 ( 1 .J.K). 83 1  I . J.K ) .84 ( I ,J.K 1,851  I , J 
2K 1.86(1. J.K). 87(1. J.K 1.88(1. J.K). 89(1, J.K! 

101  CONTINUE 
RETURN 
ENO 

6IBFTC  0ECK2  NOOECK 

SUBROUTINE  NRCELL (M.RC. PC* ZC. AL. NX I .xm.vm.zm.aa.bb.cc ) 

NX1  IS  THE  MAXIMUM  PHI-MOOE  INOEX  PLUS  I. 

RC.RC.Z'*  DEFINE  THE  CENTER  OF  THE  NEAR  CELL 
XM.VM.ZM  ARE  THE  MATH  POINT  COORDINATES. 

CELL  DIMENSIONS  ARE  R«AA.  PHI«BB.  L»CC. 

P|  THROUGH  P16  ARF  THF  VARIOUS  INTEGRALS. 

M  IS  THE  NUMBER  OF  SUBCELLS  A^.ONG  EACH  EDGE  OF  THE  MAIN  CELL. 
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e  M#M#M  is  THE  TOTAL,  no.  oi*  subcells.  M  Is  even 
COMPLEX  XP.P.O 

COMPLEX  Pi <3>  .P2<3>  .P3<3k.P4<3>.PS«3|.P6C3|.P?(3>.P8n>  .P9<3I 
COMPLEX  P|0l3>.Pl I  «3).Pl2<3» .Pi 3<3 » «PI4 1 3 > .PI  9 13 1 .Pl6<3 1 
COMMON  PI ,P2.P3.P4.P5.P6.P7.P8.P9.P10.P1 1 .PI  2 .PI 3.P 1 4.PI9.PI 6 
COMPLEX  XPP.XPO 
REAL  L.l-S 
00  90  N«t .NX  1 
PI IN)a<0.,0. I 
P2(N)«(0.,0. | 

O3(NI«(O..0. ) 

P4(N>«! 0..0. ) 

PBINla (0. <0. ) 

O6(N).(0..0. t 
P7|NI«<0. .0. I 
P8INI*(0..0. | 

P9(N)a(0..0. ) 

P|0fN|«<0..0. ) 

Oil <N|a<0..0. » 

PI2IN>*<0..0, ) 

PI3tN1«<0..0. 1 
P14(N)«(0. .0. ) 

P19<N)*«0..0. ) 

O16«N»«»0..0. ) 

90  CONTINUE 
T  l  «M 
T2-M/2 
A»AA/'T1 
B»BB/T1 
CaCC/Tl 

RAOS»«C-«T2-0.900l*A 
PH!S«PC-IT2-0.900 1*8 
l S*- IT2-0.907 ) *C 
00  100  l«I.M 
Tl-I-I 

PHl«PHIS+Tl*e 
Cl»COS(PHI 1 
SI*SIN<PM1 ) 

OO  100  J*1 .M 
Tl-J-J 

rao«raos+ti*a 
00  100  tc«l.M 
f «»«-! 
t  *LS*T1*C 

XaXM— lRAD*C05 I AL  > -L*S INI AL » )*C0SlPMI ) 

V*7M-«RA0*C0S(AL)-L*S|N( AL  >  »«SIN»PMI  ) 

2«2M-lZr*L*COS( ALI I 
RaSORT <  X*X*T* Y*Z*Z  > 

R3*P*#3 

RS*R«*9 

AR*«.28JI9*R 

PR*<-1 »*AR*AR|/R3 

P|a-ARSR3 

PaCMPLXlPR.Pl  I 

OR* I3.'AR*AR|/R9 

0I*3.*AR/R9 

QaCMPLXIOR.OI I 

AP«RAO*<*OS(ALI 

OVaAP*B*A*C/ft .2831 9 

XP*CMPLX I COS  I AR I . -S I N ( AR  >) 

XPPaXPPP 
XPO«XP*Q 
OO  100  Nat.NXl 
re«N-i 
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CN*COS ( T2+PM I ) 

SN*SINIT2*PHI 1 

CND«CN*DV 

SND»SN*DV 

PI  IN  1«P1  IN1*C1*CND*XPP 
P2IN1«P2IN1+S1*SN0*XPP 
P3IN1 »Pj  IN >4-X*X*C I  *CND* XPO 
P4INJ»P4INM-X*X*S1*SN0*XP0 
PS ( N ) -PS I N I ♦X*V*S I *CNO*XPO 
P«  IN  »  -P8IN  )*X*y*C  I  »SND* XPO 
P7 IN  » »P7 I N  >4X*Z*CNO*XPO 
08 1 N  } -P8 1 N ) +V*Y*S1 *SNO*XPQ 
09 I N | -P9 I N I ♦ Y#Y*C 1 *CNO*XPQ 
O|0(N|*P|  niNJ-»Y*Z*SND*XPO 
PJ 1 INl-Pl 1 IN)*CNO*XPP 
Pl2INl-Pl2INH-Z*X*Cl*CN0*XP0 
PI3IN»»P| 3IN>*Z*X*S1 *SNO*XPO 
O|4IN)«P14IN)*Z*Y*51*CN0*XP0 
PlSIN»«P|SIN)+Z*Y*CI*SNO*XPO 
P 1 « I N I »P 1 6 « N »  +Z*Z*SNO*XPO 
100  CONTINUE 
RETURN 
END 

•IBTTC  DECK 3  NO DECK 

^UNCTION  FHZ1 

eu-Z 

RETURN 

ENO 

804T* 

I  2  0.450  0.100 

11  2  4  6  8  10 


Fig.  26—  (continued). 
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